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Abstract. We study the directed polymer (DP) of length t in a random potential in 
dimension 1 + 1 in the continuum limit, with one end fixed and one end free. This maps 
onto the Kardar-Parisi-Zhang growth equation in time t, with flat initial conditions. 
We use the Bethe Ansatz solution for the replicated problem which is an attractive 
bosonic model. The problem is more difficult than the previous solution of the fixed 
endpoint problem as it requires regularization of the spatial integrals over the Bethe 
eigenfunctions. We use either a large fixed system length or a small finite slope KPZ 
initial conditions (wedge). The latter allows to take properly into account non-trivial 
contributions, which appear as deformed strings in the former. By considering a half- 
space model in a proper limit we obtain an expression for the generating function of all 
positive integer moments Z n of the directed polymer partition function. We obtain the 
generating function of the moments of the DP partition sum as a Fredholm Pfaffian. 
At large time, this Fredholm Pfaffian, valid for all time t, exhibits convergence of the 
free energy (i.e. KPZ height) distribution to the GOE Tracy Widom distribution 
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1. Introduction 
1.1. Overview 

The continuum Kardar-Parisi- Zhang (KPZ) equation describes the non-equilibrium 
growth in time t of an interface of height h(x,t) in the presence of noise [1]. It defines 
a universality class believed to encompass numerous models and physical systems [2]. 
The KPZ problem maps to the equilibrium statistical mechanics of a directed polymer 
(DP) in a random potential, the simplest example of a glass [3] with applications to 
vortex lines [3], domain walls [5J, biophysics [6] and Burgers turbulence [7J. Via the 
exact Cole-Hopf transformation, the DP polymer with two fixed endpoints maps onto 
the KPZ growth starting from a "droplet" initial condition, while the DP with one 
fixed and one free endpoint maps to KPZ growth from a flat initial interface. The third 
important case is the stationary growth in a finite geometry, equivalent to a DP on a 
cylinder. In the continuum, the DP and KPZ problems completely identify, and the 
free energy of the polymer is simply proportional to the KPZ height while its length t 
is proportional to the KPZ time. 

An important property of the growing KPZ interface is that it becomes, at large t, 
statistically self-affine with universal scaling exponents. In d — 1, its width is predicted 
to grow as Sh ~ t 1 ^ 3 |8J, as observed in experiments [HI EH FTT] . Furthermore, in a 
remarkable recent experiment on liquid crystals [12] it was possible to measure with 
great accuracy, not only the scaling exponents, but also the full probability distribution 
of the height field for various initial conditions. It was observed that this distribution 
depends on whether the interface starts from a droplet or from a flat configuration. 

Since the KPZ equation has resisted direct analytical solutions, progress came 
indirectly, i.e. from solving models believed to be in the same universality class. 
Analytical progress on the KPZ universality class in d = 1 came from exact solutions 
of a lattice DP model at zero temperature [13J, discrete growth models such as the 
PNG model [HI US], asymmetric exclusion models [16J and vicious walkers [T7J. An 
analog of the height field h(x, t) was identified and, in the large size limit, its one-point 
(scaled) probability distribution was shown to equal the (scaled) distribution of the 
smallest eigenvalue of a random matrix drawn from the famous Gaussian ensembles, 
the so-called Tracy Widom (TW) distribution [18], which also appears in many other 
contexts pj)]. It was found [HI [20] that depending on the boundary condition, e.g. 
droplet versus flat initial condition for the interface, one gets either the TW distribution 
i*2(s) of the Gaussian unitary ensemble (GUE) or Fi(s) of the Gaussian orthogonal 
ensemble (GOE) for droplet and flat initial conditions respectively. The corresponding 
many point distributions were identified as determinantal space-time processes, the 
Airy process, Ai 2 for droplet, and Ai\ for flat which are both naturally expressed as 
Fredholm determinants [2T] . 
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These advances gave valuable, but only indirect information on the continuum 
KPZ equation, i.e. a conjecture for its infinite t limit (termed the KPZ renormalization 
fixed point [22J ) , but no information about the approach to this limit, that is of great 
importance for experiments [H [101 HH [12] , which always probe a large but finite time 
window. This long time approach to the asymptotic value is particularly relevant in the 
recent experiments p2] which seem to confirm that the height fluctuations are indeed 
described by the GUE and GOE TW distributions at large time depending on the type 
of initial condition. 

Recently, we found [231 121] an exact solution for the continuum directed polymer 
with both endpoints fixed, i.e. the KPZ equation with droplet initial conditions, by 
exploiting the exact Bethe ansatz solvability of the replicated problem (initiated in 
Refs. [213 [2H1 I27J)- A similar result for the droplet initial condition was derived within 
different approaches [281 ES]. This solution shows that the generating function (called 
g(s) below) of the free-energy distribution is a Fredholm determinant, not only in the 
infinite time limit (the TW distribution is naturally written as a Fredholm determinant), 
but also for finite time. It thus provides an exact solution for the universal crossover 
in time in the continuum KPZ equation. It was also shown [23] that this universal 
distribution (depending on a single parameter t) describes, in the DP framework, the 
high temperature regime [23] that has remarkable universal features [30J . In the growth 
problem, this corresponds to a universal large diffusivity-weak noise limit, at fixed 
correlation length of the noise. 

Several authors then exploited the powerful replica Bethe Ansatz approach to the 
KPZ equation and the directed polymer, and a number of other quantities and initial 
conditions have been considered [311 E21 ESI EH ESI ESI EH EH1 EH]- In particular, very 
recently, the two remaining important classes of initial conditions have been solved using 
replica, the flat initial condition [35] and the stationary class [37]. Other remarkable 
recent results concern the distribution of the position of the DP endpoint (101 E] and 
have been obtained by completely different methods. Some recent more mathematical 
achievements have been reported in Refs. [12] (see also Ref. [13] for a review). 

1.2. Quantities of interest, strategy of the calculation and main results 

The aim of this paper is to provide a detailed account of our recent solution of the KPZ 
equation with flat initial condition [35], i.e. the directed polymer with one free end and 
one fixed. The continuum model has natural space and time units and we define 

A = \(cH/T^, (1) 

the single dimensionless parameter which describes the crossover from the diffusive 
regime at small time t to the glass fixed point regime at large t. Here T is the 
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temperature, c the disorder strength and t the length of the DP. In terms of KPZ 
parameters this is T = 2v and c = D\q. All parameters are defined below. Note that 
([!]) is exactly the same parameter A as in our previous work [23] (and in [21]) on the 
fixed endpoints problem. 

Our aim here is to characterize the probability distribution (PDF) of the free energy 
of the DP 

F = —T\nZ = TXf , (2) 
equivalently, of the height field of the KPZ equation 

^h = ]nZ = v t + \£ t . (3) 

To achieve this, we first calculate by means of the Bethe ansatz the disorder average 
(denoted here by overbars) of the integer powers of partition function Z n as function 
of time. We then avoid the explicit n — y limit by introducing the generating function 

gx ( s ) = e-*- XsZ = 1 + V 1 — ' Z» = exp(-e- A ( s +/)) , (4) 

/ — ' n! 

n=l 

where the subscript A is often omitted below. Once g\{s) is known, the PDF of the 
rescaled free energy P(f) at large time (i.e. A — » oo) is immediately extracted as 

+oo 



lim g x (s) = 6(f + s) = Prob(/ > -s) = / dfP(f) . (5) 

A^oo J_ s 

and it can also be extracted at finite time (e.g. via the procedure described in [23]). 
In this paper 9{x) denotes the Heaviside function. The same generating function was 
introduced and computed in [23] for the fixed endpoint problem when it takes the form 
of a Fredholm determinant (FD) at any time t (any A). Here the calculation is much 
more involved but we finally find that g\(s) takes the form of a Fredholm Pfaffian (i.e. 
the square root of a Fredholm determinant) for any A. We then show that at large time 
this Fredholm Pfaffian, giving the free energy (i.e. KPZ height) distribution, converges 
to the GOE Tracy Widom distribution, i.e. the scaled random variable £ t = — / defined 
in (]3j obeys 

9oo ( S ) = lim Prob(6 <s)= F 1 {s) . (6) 

t— >oo 

As we will see, the calculations are far more complicated than for the DP with fixed 
endpoint. The first step, i.e. obtaining explicit expressions for the integer moments Z n 
is much more involved. In the fixed endpoint calculation one could directly consider 
the infinite space L = +oo and perform a summation over the string states, which 
have an easy structure in that limit. Now however we need to calculate the spatial 
integrals over the Bethe eigenfunctions, and these are not nicely convergent on the 
infinite space (one easily gets either zero or infinite result!). The first approach we try, 
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which we call direct approach is to consider (i) a periodic system of finite size L, for 
large L. The correct treatment however requires the use of the Bethe equations which 
enforce periodicity. This approach quickly becomes intractable, and we have been able 
to calculate only small integer values of n (we explicitly report only up to n = 4), or to 
sum for all n but at fixed number of strings n s = 1,2. A structure in the contributions 
however starts to appear, confirmed later by other methods. The second approach is 
to consider (ii) a small finite slope w > in the KPZ initial conditions, i.e. a wedge 
initial condition. There the limit L = +oo can be taken from the start and the desired 
integrals can be obtained by considering poles in 1/w. The symmetric wedge allows to 
make some progress, in particular we identify a contribution which we missed at first 
in the direct approach, which can be identified as a deformed string in the large L limit. 
Hence, both methods taken together give reliable results for small values of n and/or 
n s . The true progress however comes from considering the half-space wedge. By taking 
a combination of limits (a) large distance from the wedge (b) small w, we are able to 
compute systematically all contributing poles in 1/w and derive their Pfaffian structure. 
The rest of the calculation relies on further tricks to perform the summation over all 
string states in the form of a Fredholm Pfaffian (FP). This FP takes various forms, 
some of them can be useful for numerical evaluations, while one is particularly suited 
to take the large time limit and to show convergence to the GOE TW distribution. 



1.3. Organization of the manuscript 



The paper is organized as follows. In Sec. [2] we introduce the DP and KPZ models, the 
various boundary conditions, and the notations employed in this paper. In Sec. [3] we 
recall the needed details about replica and the Bethe ansatz solution for the attractive 
Bose gas. The starting formula for the moments and the string partition sums are 
given in Sec. 32 and 3J3 and the Airy trick which allows to perform divergent sums is 
recalled. In the following sections we start our new computations. As a first step, in Sec. 
[4J we attempt to calculate Z n directly for the flat initial condition. Although it does 
show a nice structure emerging, the actual calculation is very cumbersome and only 
results for n = 2,3,4 and n s = 1,2 are reported. The exact result for n s = 1 already 
reproduces the correct right tail of -Fi(s) (i.e. left tail for P(f)), as discussed in Sec. 



4.2.1 In Sec. p5lwe move to the wedge initial condition. The full-space or symmetric 



wedge is discussed only in the Appendix D but as we show it allows to reproduce flat 
or narrow initial conditions, respectively. We thus proceed with the half-space wedge 
and analyze its pole structure in Sec. 573] by taking proper limits. This leads to the 
full analytical formula for Z(n s , s) which is presented in Section 5.5, and checked to 
agree with all previous small n, n s results by the two other methods, as well as with 
a non-trivial sum rule at t = 0. The reader uninterested in the technicalities involved 
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in deriving the expression of Z(n s , s) can skip sections [4] and Appendix D In Sec. [6] 
we proceed with the Airy trick to derive the final expression for g(s) as a Fredholm 
Pfaffian. It involves on the way some new Pfaffian identities which allow to perform 
the summations over the number of particles in the strings. The explicit expressions 
for the kernels are given in Sections 6.3 to 6.5| The large time limit is then studied in 
Section [7] and the convergence to the GOE Tracy Widom distribution is proved. Section 
[8] contains concluding remarks. In a series of appendices we report the most technical 
parts of the calculations. 



2. Models and boundary conditions 

2.1. Directed polymer in the continuum 

Here we study the continuum model for the directed polymer in d — 1 + 1 dimension in 
a gaussian random potential, defined by the partition sum Z(x, t\y, 0) of paths x(t) e R 
starting at x(0) = y and ending at x{t) = x 

PX [t^j — X 

Z(x,t\y,0)= / Dxe^ti d ^) 2 +V(^M], (7) 

Jx(0)=y 

which is the solution of the Feynman-Kac equation 

d t Z = -H d Z = T -V\Z - ±V(x, t)Z , (8) 

with initial condition Z(x, t = 0\y, 0) = 5(x — y). In the (^-correlated continuum model, 
the random potential V(x,t) is furthermore chosen centered Gaussian with correlator 

V(x, t)V(x', t) = c5(t - t')5(x - x). (9) 

By a change of units, i.e. by a rescaling x — > T 3 x and t — > 2T 5 t one can eliminate the 
temperature T: below we work in these units hence set T = 1 everywhere. In addition 
it is also possible to also set c = 1 by a further change of units x — > x/c and t — >• t/c 2 , 
a freedom that we sometimes use below (when explicitly indicated). 

As discussed in details in [30j[23], the 5-correlated continuum model ^ describes 
the universal high T limit of (i) the DP on a lattice (ii) a continuum DP model, but 
where the disorder has a finite correlation length rj along x (by contrast r/ = in Q). 
In both cases the description by the 5-correlated model becomes exact when the large 
scale which appears at high temperature is much larger than the small scale cutoff of 
the model, e.g. T 3 /c ^ rj in case (ii). 
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2.2. KPZ equation in the continuum 

An important motivation to study the DP problem comes from its relation to the KPZ 
growth equation for the height field h(x, t) 

d t h = isV 2 h + ^\ (Vh) 2 + Z(x,t), (10) 

in presence of a noise t), where v is diffusivity and the strength Ao of the non-linear 
term has no relation to the parameter A. The Cole-Hopf mapping 

Z(x,t) =e^ x ' t \ T = 2v, V(x,t) = -\ £(x,t), (11) 



maps (|8j) into (|10j). More precisely, the exact relation in terms of the initial condition 
h(x, t = 0) for the KPZ interface is 

e £h(x,t) = J dyZ(x,t\y,0)e^ h{y ' t=0) . (12) 

where Z(x, t\y, 0) is the DP partition function given by (J7|). At this stage, these relations 
are valid for arbitrary V(x, t) and £(x, t) (and also in any dimension but here we study 
d = 1). In the standard KPZ problem the noise is also assumed to be a centered 
Gaussian white noise with 



£(x, t)£(x', t>) = D5(x - x')5{t - t') , (13) 

and so it corresponds to the 5-correlated continuum DP (|9]), upon the further 
identification c = -DAq. We mwntion that there are some mathematical difficulties in 
defining KPZ equation with white noise [2HJII2], but that they are usually surmounted 
by defining it precisely as the Cole-Hopf transform of the DP problem, which does not 
suffer from such problems. 

By correspondence with the DP (see above), the white noise KPZ equation 
describes a universal limit of growth models. For instance, if we consider a continuum 
KPZ equation where the noise has finite correlation lengths (along x and t) it will be 
described by the white noise limit whenever the characteristic time t* = 2(2u) 5 /D 2 Xq 
and space x* = \Jvt* scales are much larger than these correlation lengths [fj Hence it 
corresponds to the high diffusivity v or weak noise D limit. Similar statements can be 
made for discrete growth models although we will not elaborate here. 

2.3. Boundary conditions 

It is well known that the critical exponents are independent of the chosen initial 
condition [13] and so they have an high degree of universality. By contrast, the 
distribution of the fluctuations of the height field h(x, t) (i.e. the probability density 

| in that case D is defined as the space time integral of the two point correlator, see Refs. [3D1 [23] for 
the same statement on the DP 
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distribution of the free energy of the directed polymer) depends on the initial conditions 
h(x,t = 0) (or at least some features of them) even for large time and so the 
initial condition specifies different universality classes (see [13] p. 16 for a detailed 
classification). Three main categories of initial conditions have been identified that 
are (i) the sharp wedge (also called droplet) initial condition, (ii) the flat one and (iii) 
the random one (which includes the stationary one). For large time, the probability 
distribution functions are expected to crossover to distinct universal distributions, some 
related to the Tracy- Widom distributions for the largest eigenvalues of various random 
matrix ensembles. In addition to the infinite time behaviour, for the continuum KPZ 
equation with white noise (as defined above) each main class exhibits also a one 
parameter family of universal finite time crossover distributions (indexed by A) and 
here we will mainly study the one associated to the flat initial condition. 

A slightly more general, and particularly useful initial condition that we also study 
here is the so-called wedge initial condition 

^h wedge (x,t = 0) = -w\x\, (14) 

that has the important property to interpolate (for any t) between the infinitely narrow 
wedge for w — > oo (sharp initial condition) and the flat interface initial condition for 
w — > + , thus having as limiting cases two of the three important random matrix 
ensembles. In addition, as discussed below, for fixed w but infinite t it converges to the 
fixed endpoints (droplet) class, as physically reasonable. 

In our previous work [23] we studied the sharp wedge initial condition for the KPZ 
interface, that corresponds to the DP with the two ends fixed (say at xo = 0), i.e. 
Z = Z(0,t\0, 0), and obtained the distribution of InZ for any time. In this paper we 
are interested in the flat interface initial conditions of the KPZ growth problem, i.e. 
h(x, t = 0) = 0. This corresponds to a DP where one end has been fixed and the other 
is instead free, so that the partition function of interest is 

Zfiat(M) = J dyZ(x,t\y,0), (15) 

i.e. the sum of all directed paths which join a given point x(t) = x to any point y 
on the line at t — 0. As already stated, in terms of the wedge initial condition, this 
corresponds to the opposite limit w — > + . We will also discuss here the general case, 
i.e. for any w > 0, which interpolates between the two problems, but we are still unable 
to provide a full solution. 

In order to achieve the solution of the flat initial condition, we first attempt to 



solve directly the case w = 0. Although the integral (15) is perfectly convergent, due to 
the fast decay of the diffusion kernel (stretching energy of the DP), divergences appear 
in the Fourier representation of the Bethe states in the infinite space problem, and 
force us to start from a finite (periodic) system of length L along y, and take the limit 
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L — > oo only after the calculations of the moments of the partition function Z(x, t) 



in (15). This way of proceeding is plagued by further technical difficulties that allow 



us to get only partial results. For this reason we move to the full-space wedge initial 



condition ( 14 ) with DP partition sum 

/oo 
dye- w ^Z(x,t\y,0), (16) 
-oo 

with Z(x,t = 0) = e~ w l x l. In this case all divergences are canceled by the exponential 
decay of the initial condition and it is possible to work directly in the thermodynamic 
limit L = +oo. We reproduce all the results found in the previous approach and extend 
them. However, some other technical problems make difficult a full solution. While it 
does not seem impossible to overcome these technical issues, we have found much easier 
to consider the wedge initial condition for (left) half-space problem with partition sum 

Z hs , w (x,t)= f° dye w yZ(x,t\y,0), (17) 



with Z(x, t = 0) = 9{— x)e wx . As the physical intuition suggests, this half-space problem 
has the important property that at fixed w = + it interpolates between (i) the narrow 
wedge initial condition for x — > +oo, since the polymer is then stretched and (ii) the 



flat, full-space initial condition (15) for x — > — oo, because very far from the origin 



the polymer does not feel the presence of a boundary. It is easy to show this from 



the statistical tilt symmetry (STS) of the problem (see Appendix A). This implies the 



relation for the moments of the partition function for the half-space problems 



Z»(x,t) = e"*-ZZ + j s (0,t). (18) 

Thus changing the endpoint is the same as changing w, up to a simple additive piece in 
the free energy. The problem with w = and x > is the same as w = t| and endpoint 



at 0, and in the limit x — > +oo it thus identify with the sharp edge problem. Eq. (18) 



also implies that Z™(x, t) = e nw ( x + wt ) Zq(x + 2wt, t). Hence if x — > — oo faster than 2wt 
tends to +oo the problem becomes a full space one with an edge far away, as claimed 
above. The only difference is a uniform shift 2wt in space and w(x + wt) in the average 
free energy. 

Thus the final solution for the flat initial condition will be achieved by considering 
the limit 

lim Z hs (x,t) = Z flat (x,t). (19) 

x— >— oo,w— >0 

The limit being taken as discussed above. Note that this limit still has x dependence 
in a given disorder realization. 
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3. Quantum mechanics 

The calculation of the n-th integer moment of a DP partition sum can be expressed as 
a quantum mechanical problem (see e.g. [251 EZ])- Upon replication of (j7l), averaging 
over disorder and using the Feynman-Kac formula, one finds that the moments 



Z n := Z(xit\y 1} 0)..Z(x n t\y n , 0) , (20) 

satisfy 

dtZ n = —H n Z n , (21) 

where H n is the (attractive) Lieb-Liniger (LL) Hamiltonian jH] for n particles with 
interaction parameter c = — c 

^ = -E^ + 2c E (22) 

j=l X<i<j<n 

In other words the Z n can be written as (imaginary time) quantum mechanical 
expectations: 

Z n = ( Xl , ..x n \e- tHn \ yi , ..y n ) . (23) 

At this stage, and to compute Z n for arbitrary endpoints one would need all eigenstates 
of H n , not only the symmetric ones (bosons). As explained below however, for the 
observables computed here we will need only the bosonic ones. 

This many body quantum mechanical model is solvable by means of Bethe ansatz 
[44J. This consists in an educated guess for the symmetric eigenstates of the many- 
body wave-functions ^^(xi, ..x n ) = (xi, ..x n \fi) that are superpositions of plane waves 
jB] of the form 

»„(*„. .*„)=5>n a p = n i i - ic T ( '-x x ° ) )- (24) 

P j=l n>p>a>l P P Pa 

The sum runs over all n\ permutations P of the rapidities Xj. As written, these Bethe 
wave functions are not normalized and satisfy 

(x- ■ • x\(jl) = ty^x, ..x) = n\e ix ^ Xa . (25) 

The rapidities must be all different, i.e. X a ^ \p for a ^ (3 (the wavefunctions are 
antisymmetric in exchange of 2 rapidities, hence vanish if 2 are equal). When working 
with periodic boundary condition, also ^ ^(xi, ..x n ) must be periodic in any of the 
variables and this forces the set of rapidities {A} to be solution to the equations 



e^ L = TT ^~^ + * C , (26) 
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that are the Bethe equations for the Lieb-Liniger model. In particular these imply 
e'^ lA ° = 1 and hence the total moment K = ^" =a A a = 2irp/L is quantized in 
terms of integer numbers as it must be. 

The solution for the attractive case has been described soon after the Lieb-Liniger 
repulsive solution by McGuire [HJ. We report here the main ingredients to understand 
the structure of the wave functions. For c = — c > the Bethe equations are 
i\ a L _ TT A a - - ic 

Consider now a complex rapidity X a = A + if]. The Bethe equation for this rapidity is 

e *X a L = e iXL-r,L = TT Att-Afl-ic 

f} X a -Xp + ic 

We consider n to be finite and L — > oo. If rj > 0, we have e~' nL — > on the left-hand 
side. Looking at the finite product on the right-hand side, we conclude that there must 
thus be a rapidity X a > such that X a i = X a —ic+0(e~ vL ). On the other hand, if 77 < 0, we 
have e~ vL — > 00 on the left-hand side, and there must thus be a rapidity X a > such that 
X a > = X a + ic + 0(e~' ?? ' L ). Furthermore if A = X a is a solution also A* (if A 7^ A*) must 
be for some (3 7^ a as evident by taking the complex conjugate of the Bethe equations 

Thus the rapidities solution to the Bethe equations for c < and L — > 00 arrange 
in object formed by an arbitrary number of particles m < n with the same real part 
and imaginary parts that differ by integer multiples of ic and that are symmetric with 
respect to the real axis. These objects are called strings. A general eigenstate is built 
by partitioning the n particles into a set of n s < n strings formed by rrij > 1 particles 
with n = Y^j=i m r The rapidities associated to these states are written as 

%c 

X^ a = kj + —{mj + l- 2a) + i5 J ' a . (29) 

Here, a — 1, ...,rrij labels the rapidities within the string j — 1, . . .n s . 5^ a are deviations 
which fall off exponentially with system size L. As it should be clear, perfect strings 
(i.e. with 5 = 0) are exact eigenstates in the limit L — > 00 for arbitrary n. The 
(unnormalized) wave function for a string of m particles (which in the following will be 
simply called an m-string) is 

..x m ) = nlexp ^2 \ x i- x j\)> ( 30 ) 

l<i<j<m 

and so each string represents a bound states of rrij particles with its own energy and 
momentum given by 

c 2 

Kj = rrijkj, E, = rrijkj - — m j( m2 j ~ ( 31 ) 
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The fact that all deviations are exponentially small in L implies that we can obtain a 
much simpler set of equations involving only the string momenta Kj reported in Ref. 
|46j, that in the limit of L = oo gives vanishing coupling between strings. Thus for all 
our aims, the only effect of the (attractive) interaction is to produce a set of bound 
states. When instead interested in the directed polymer in a finite geometry, other 
methods to tackle exactly the Bethe equations should be used, as done e.g. in Ref. 

The general eigenstate consisting of partitioning the n particles into a set of n s 
strings formed by rrij > 1 particles with n = Y^jLi m j have momentum and energy 

i=i j=i ^ ' 

Clearly, by minimization of the energy, the ground-state is the single string built with n 
particles and having zero total momentum, while excited states are obtained by either 
giving momentum to the ground state, or by partitioning it into smaller strings to which 
individual momenta can be given. 



3.1. The replica method for the flat and wedge initial condition 

In this paper we are interested in the moments of the partition function for the flat 
initial condition in Eq. ( 15 ) . By definition in the replica approach these are 

(yi 



J flat 



[x,t)< 



3=1 



Vn\e 



-tH n 



\x . . . x) . 



(33) 



As compared to (23) we have, for future convenience, exchanged the starting and the 
end points of the polymer, using that Zfl at (x,t) n is real and H n hermitian. Assuming 
completeness of the string states, we can insert a resolution of the identity over the 
Bethe states (representing all the possible allowed rapidities) and obtain 

(yi ■ ..y n \(i){fi\x...x) 



J flat 



x,ty 



V ^fo"*) r ~tE u 

^ IMI 2 



-tE u 



/i • • • y n ) 



(34) 



Notice that it has been possible to use the solution of the bosonic problem only because 
the state \x . . . x) is symmetric for exchange of particles (replicas). Thus it has vanishing 
overlap with anti-symmetric states representing fermionic wave functions which are not 
included in the above Bethe ansatz solution. 



In Eq. (34), the energy is given in Eq. (32) and the equal points wave function 
^*u{x, ..x) in Eq. (25). The norms of the Bethe states within our normalization have 
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been calculated in Ref. jl6] by adapting the famous Gaudin-Korepin formula [50J for 
norms of Bethe states to the case of strings. The final result can be written as [46J 



1 c" x „ lT 4 1 ^r, tt 4(fcj - A;,-) 2 + (m; - m i ) 2 c z 

- ,. r , $ fc,m I J — 3 , $&,m = ] — ^. 35 



It is then evident that the only ingredient missing in Eq. (34) is the integral of Bethe 
wave function for string states that will be one of the main tasks of the following. 

It is straightforward to write down the corresponding formulas for wedge initial 
condition both in full- and half-space problem. We have simply 
\d*(x t) r w 



where we defined 



/w poo . ,v . 

^ := / ( JJrf%)e" w ^ n = ll2/ill I r At (yi, ..y n ) , for the full-space problem, (37) 
J — OO j = l 

/w pQ . n . 

^ : = / ( J\dyjje w '^*= lVi '^ li (y 1 , ..y n ) , for the half-space problem. (38) 

Note that in the limit w — > +oo one should recover the DP with both endpoints fixed 
(one at x the other at 0) i.e. the droplet initial condition, by the simple replacement 
e ~ w \y\ _). —5(y) for the full space and e wy 6(—y) — > —6(y) for the half-space model. The 
factors 2/w and 1/w lead to an additive normalization ~ Inw in the free energy. 

3.2. Starting formula for moments 

To summarize we have the starting expression for the moments 

z(x~tr = y ° y f\— 2 

v ' ^ n s \(cL) n ° ^ Urn 2 

1s=l (mi,...m„J„ 3=1 J 

x Z)( / m ] fl e lra '" ra ^" mjfc ' i "™ jfc) , (39) 

where J ^ stands for the various cases, and is a function of [k, m\. Here (mi, . . . m n J n 
stands for all the partitioning of n such that X]j=i m j = n with m j > 1- I n the limit 
L — > oo the momenta sums become continuous and one can use that the string momenta 
mjkj correspond to free particles i.e. ^2 k . — > rrijL J We will however also often 



make use of the discrete sum (39) below whenever the momenta are constrained to 
special values, e.g. to vanish. 
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3.3. Starting formula for the generating function 
Let us recall that the generating function 



7(t t\ n 1 
g(s) =1 + J2 ^£L(-ire-^ = 1 + J2 ±Z{n a , s) , (40) 

n» = l 



admits an expansion in the number of strings 
Z{n„s)= (cLWV.m,-)!!!^ 

m 1 ,...m„ a =l K ' KL ^J J J j=l 3 

r Us 2 

x / *^[^]n e(raj3 " mj) ^"" ji)2i " imjfcrAmjs ' ( 4i ) 

obtained by permuting the sum over the number of particles and strings, which allows 
to free the constraint on the summations over the number of particles, as in a grand 
canonical ensemble where As acts as a chemical potential. 



As discussed in [231 [24], while the formula (39) is perfectly well defined, because of 
the exponential cubic divergence of the series, the formula (41) should be taken in some 
analytical continuation sense, i.e. it is valid as a formal series in t. This continuation 
is performed in the following using the famous Airy trick which for Re(w) > ensures 

dyAi(y)e yw = e w3/3 . (42) 

the summations can then be later carried at fixed y. 

4. Direct approach to the flat initial condition 



From Eq. (34) we need to evaluate the integral 



n^iKfa ••■!/»), ( 43 ) 

over all possible string wave functions. In order to do so let us start by considering the 
auxiliary integral over the principal domain PD = {{yj}\0 < y\ < y 2 < ■■ < y n < L} 

B n (A 1 ,..A„) = / d yi ..dy n e lX ^ + - +lX ^ = i n { V(-1)V l ^=»-j+i*< J- + J-), (44) 
Jpd CLjbj a Q b 

where we introduced 

j n—p+l re— 1 n—j 

p=l i=n—j+l p=j i=n—p 
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Clearly this formula is true if and only if the set {Ay| i<i<j< n } with 



p=i 



does not contain any zero, else the corresponding term in the sum diverges. Using the 



symmetry under permutation of the yj we have for a general Bethe state (24) 



[d n y^,(y) = n\J2B n (X Pl ,..X Pn ) ]J (1 - — ). (47) 



n>^>fc>l 



Fixing the values of Xj by using the Bethe equations (26), it is possible (but 
cumbersome) to show that this combination always vanishes: 

dVMy) = 0. (48) 

This result shows that in order to have a non- vanishing contribution of the Bethe 
wave function, we need to require that at least one of the A^ should be zero (up 
to permutation). To obtain the result in that case one needs to (i) enumerate all 
possibilities for a subset of Ay to vanish (ii) perform limits in the above expression. 
Below we give some partial results, with few details, being assumed that in each reported 
case we have searched all such possibilities and performed the necessary limits (which 
is tedious). 

In order to search for the general structure of the states that contribute non- 

trivially to this problem, we first consider the calculation at fixed n, i.e. the average of 

the integer moments of the partition function of the directed polymer. Next, we also 

obtain some results for fixed number of strings n s = 1, 2. For actual calculations a 

useful formula is: 

r°° i n i n 
\ dLe- sL B n (X u ..X n ) = -T\ , /x i = y)A J , (49) 

JO S ■ i S l t L i ■ ■ 

Upon summing over permutations, Laplace inversion and using the Bethe equations 
(BE) we obtain the results discussed below. One can see on this formula that multiple 
poles in 1/s lead to powers of L. 

In the rest of this section, to have a light notation, we remove the subscript "flat" 
from the partition function and we fix (without lost of generalities) x = 0. Thus 



everywhere in this section Z n stands for Z^ t (x,t) n in Eq. (34). 



4-1- Fixed n results 

We show how to obtain from the previous expression Z n for n = 2, 3, 4. An important 
check that will be a guide to check the correctness of our manipulations and to 
understand if we considered all the possible Bethe states is that for t = all the results 
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should just give the trivial expectation value of the flat condition h(x,t = 0) = 0, i.e. 

= 0) n = I- 

4-1.1. n = 2. This is the simplest possible case. Clearly there are no Bethe states 
with Xj = 0, thus the only possibility is Ai = — A2. For these states, simple algebra 
leads to 

d 2 yV(y) = -2cL^ . (50) 

There are two possible string states satisfying Ai = — A2, i.e. 

• one string of two particles (2-string state) with zero momentum (n s = 1); 

• states with two particles (1-string) with opposite real momentum Ai = — A2 = k 
(n s = 2). 

We can now write the contribution of each of these states to Z 2 using the general 



formula (39). Summing up these contributions and integrating over k we obtain the 
exact result for the second moment 

^ = + I ^(-^4)7^-^ (51) 

2\(cL) 2 J 2tt k\ Akj + c 2 

= 2e& - 4c f ^^— 2 e- 2 ^ = e 2X \l + erf(A 3 / 2 v / 2)) , (52) 
with erf (2) = 7r~ 1//2 J z dte~ t2 and A 3 = c 2 t/A. Notice the 1/2! factor in front of the 



integral comes from the 1/nJ in (39) and avoids double counting of the states. For 



t — (i.e. A = 0) we recover Z 2 = 1 as a non-trivial check of the calculation. 

4-1.2. n = 3. Only a subset of the states with some Ay = provide a non- vanishing 
contribution to Z 3 . We report here only the calculation for the states with a non-zero 
contribution that are 

• A 3-string state of zero momentum (n s = 1) for which Xf = X 2 = —c 2 and A2 = 0. 
The integral is 



36L 



d*y*(y) = . (53) 

One 2-string with k = 0, i.e. A 2 . = — c 2 /4 and one particle with k = (n s = 2), 
with integral 

J d 3 y^(y) = 72^. (54) 

Three 1-strings with (Ai, A2, A3) = (kx, k 2 , k%) = (0, k, —k) (n s = 3), with integral 
j3 T . . 6cL 2 (c 2 + A; 2 ) 

d 3 yV(y) = y — '- . (55) 
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We now sum up these three contributions to Z z using the general formula (39) and 
displaying all factors for clarity 

~Z 3 = Ae 2cH + — 2-(-72—) ° 2//4 e^ cH (56) 

+ 2\(cL) 2 A [ ' c j 9c 2 /4 e W 

o c 3 T f dk6cL 2 (c 2 + k 2 ) . k 2 4k 2 2k 2, 
+ 3 ,,_ rsn L / y — -(— -) 2 — -e~ 2k 1 



3!(cL) 3 J 2tt k 4 y k 2 + c 2J Ak 2 + c 2 

a 2c 2 t o ic 2 t o f dk 4ck 2 _ 2k 2 t 
4e 2c * - 2e2 c - 3 / — — ; rr , , ^e 2k 1 



2vr (k 2 + c 2 )(Ak 2 + c 2 ) 

where the factor 3 in front of the integral in the last line is the symmetry factor for the 
number of distinct triplets (0, k, —k). Integrating over k we obtain the exact formula 
for the third moment: 

= 4e 8A3 - 2e 2A3 - 2e 8A3 erfc(A 3 / 2 2 v / 2) + e 2A3 erfc(A 3 / 2 v / 2)) , (57) 

with erfc(x) = 1 — erf(x). For t = it reproduces correctly Z 3 = 1 

4-1.3. n = 4. Let us consider now n = 4 and report only the states giving a non-zero 
contribution. We will give only the final result while the detailed contributions and the 
integrals are given in |Appendix B The contributing states are: 



n s = 4, four 1-strings with k% = — ki and k^ = —k%. 

n s = 3, two 1-strings of opposite momenta k and one 2-string of zero momentum 
(i.e. total state — i|, k, — k)). 

n s = 2 with two 2-strings, of opposite momenta, i.e. (k + ic/2,k — ic/2,—k + 
ic/2,-k -ic/2). 

n s = 2 one 3-string and one 1-string both of zero momentum, i.e. (ic, 0, — ic, 0). 
Note that from very general arguments this state should not exist, because two 
rapidities have the same value. However, this state exists and is the infinite L 
manifestation of a deformed string (see for a discussion of the topic Ref. |51j). In 
the specific case, the deformed string originating this state is (ic+5\, e\, —ic+Si, e 2 ), 
where are the string deviations discussed previously. For any finite system e x ^ e 2 
and so the state is allowed and becomes pathological only in the thermodynamic 
limit. As discussed in Ref. |HT] these deformed strings are essential to ensure the 
completeness of the string states. Here we missed this term at first, but found 
it later by considering the wedge initial condition, which allows to treat directly 



L = +oo, in the limit w = + . This is detailed in Appendix B 
• n s = 1, i.e. a single 4-string of zero momentum. 

Summing up all the non- vanishing contributions we have 

^4 = 8e 20A3 - 8e 8A3 - 4e 8A3 (-2erfc(2A 3 / 2 ) + e 12A3 erfc(4A 3 / 2 )) 
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f 2x+l _ V2x+l } e -8X 3 x 

An(Ax(x + 3) + 5) 



- 4e 18A3 erfc(3v / 2A 3 / 2 )e 2A3 + 48 f° dx^^ — 

Jo 

1 + 12^?A 3/2 + 51.6394A 3 + . . . , 



(5f 



that for A = gives 1 as it should. The contribution of the deformed string is essential 
in order to recover Z 4 = 1 at t = as a reflection of its importance for the completeness. 

4- 1-4- First conclusions. Although some rules seem to emerge from the study of 
n = 2, 3, 4 as to which are the contributing states, we will not pursue further here 
the direct approach to arbitrary n. It becomes indeed very tedious for higher n since 
the BE must be used to show that some states will not contribute (in an a priori non- 
obvious way). At this stage we found that contributing states contain either strings 
of zero momenta or pairs of strings of opposite momenta. This rule will appear much 
more clearly and easily in the wedge initial condition studied in the next Section. 
More precisely it will appear that the contributing states are obtained by splitting the 
n s = 2N + M strings into N pairs of strings of opposite momenta and M single strings 
of zero momentum. It would be nice to prove this rule in the direct approach, but we 
have not attempted it. Of course, once the rule is known it may be possible to pursue 
the calculation in the direct approach, but we leave this task for the future. 

To close the study at fixed but large L we now show that some interesting but only 
partial results can also be obtained for arbitrary n but fixed n s . 

4-2. Summation at fixed number of strings n s 

We can also perform summation over n at fixed number of strings n s and obtain the 
first terms in the generating function ( 40|41 ). 



4-2.1. One string contribution n s — 1. It should be clear from the previous section 
that for n s = 1 the only non-zero contributions come from the n-string wavefunctions 
with k = 0, i.e. the ground state, that is 

^(xi, ..x n ) = n!e"^'<j lXl ~ Xjl , (59) 

whose integral is calculated by standard methods 

cm-\ 2 T 

d n y*»(y) = —^r- ■ (60) 



Using Eq. (39), this gives the contribution of the ground state to the moments for 
arbitrary n 

Z^\ na=1 = 2 n - l e ^ {ni - n)t . (61) 
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It is convenient to perform the usual shift in the free energy i.e. define Z = Ze~ c2t ^ 12 , 
hence g(s) = g{s + y-) and to study from now on g(s). This is equivalent to set 

^ = -y, (62) 
in Eq. ([3]). Inserting the above expression in Eq. (41) we obtain the very simple result 

Z(l, s) = V ^ ' 2 m - 1 e m3 ^- Xms . (63) 

^ ml 

m=l 

Performing now the Airy trick we have 

Z(l, S) = J dyAl{y) t^l 2 m-l e Xmy-Xms = ^J dyM{y + ^-fc* _ y ^ (g4) 



m=l 

that for t — > oo is 



lim Z(l,s) = --[ dyAi(y + s) = -- I dyAi(y) . (65) 



A ^+°° 2 Jy>0 ' ^Jy>-s 

Analysis of Z(l,s) at infinite X: although quite simple, it is remarkable that Z(l,s) 
at infinite A already reproduces the first term in the asymptotic expansion of the 
GOE Tracy-Widom distribution F\(s). Indeed, as explained in more details in 
Section [7j F\(s) = Det[J — B s ] is a Fredholm determinant of the kernel B s (x,y) = 
8(x)Ai(x + y + s)6(y) and as such it can be formally expanded in powers of B s , i.e. 
Fi(s) = l — TrB s + 0(Bg). The one string contribution (i.e. the ground state for any n) 



in Eq. (65) is exactly — TrB s . It turns out that this expansion is organized exactly in 
the same way as the expansion for large negative s, which gives P(f) for large negative 
/ (deep states), i.e one finds, from the asymptotics of the Airy function at large positive 
argument: 

My) « My)*-*** . My)-^ymf o ^f^ Ai - (66) 

(fa-l)(6n-5) 
72n 

i.e. An = 5/72, Ai 2 = 385/10368, .., that (for infinite A): 

Prob(/ > -s) w Probi(/ > -s) = 1 + Z(l, s)«l| /i(-s)e"i ( - s)3/2 , (68) 

with h(y) = -^^(7^-1-1* with C = 1 and C k + |(fc - 1)^ = 
(— l) fc 74ifc(2/3)~ fc for > 1. We recall the discussion in Ref. [23J showing that one 
can construct an independent string approximation which gives a bona-fide cumulative 
probability gi n d(s) = e z ^' s ^ with the same asymptotics as the exact one for large s. 
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Analysis Z(l, s) at finite X: For any A, the moments (61) can be formally written as 

Z»\ laMng = J dfP 1 (f)e- Xfn , p 1 (f) = -Ai(-f + — ). (69) 

Clearly P\{f) is not a positive function and cannot be a probability distribution: it 
should not be indeed, since it contains only contributions of n s = 1. However, for large 
negative / we expect from the above considerations, at least at large enough A, that it 
gives the correct tail 

P{f) «/-*-«, Pi(f). (70) 

4-2.2. Contribution of two strings, n s = 2. In this case we find that there are only two 
kinds of contributions: (i) those with the two strings with mi = vri2 = m with non-zero 
opposite momenta (ii) those with mi ^ 777,2 both forced to have zero momentum. Thus 
we write Z(2, s) = Z a (2, s) + Z b (2, s) that we now calculate. 

Z a (2,s) is generated only by the states with n s = 2m even that must be divided 
in two strings of the same size m. Denoting with k the opposite momentum variable 
of the two strings (which have total momenta ±mk), we find that the integral can be 
written as 

2^(2^-1)! e 



after some heuristic manipulations reported in Appendix C We have also written the 



norm factor <3> corresponding to two equal size strings. Inserting this expression in Eq. 



(41 ) we get 



1 \ m f rlh I ,3 

1 1 -2m-l / Ufi/ TT e -2mk 2 t+2^-m' i -2m\s 



Z a (2,s)=J2—c 2m - 1 

^-^ m 2ti ±± 



m=l 

00 



27r- LJ -A; 2 + c 2 j 2 /4 
7=1 J ' 



E l n 2m f ^ TT \ -8mA 3 fc 2 +2^m 3 -2mAs frynS 

m J 27rll4P+j2 e ' " [U) 

m=l J j=l J 

where we have changed k 2 — > c 2 k 2 and used c 2 t = 4A 3 . 

For mi 7^ 7772 the calculation is sketched in the Appendix C . The final result is 



$ = f^_\ mi+m2 - 2 m 1 777 2 ( 777 1 + m 2 ) 2 r(mi + 777 2 ) 1 x lnin ( wtl , mil ) ^2 , ?3 x 
\C/ \mi — 7772|r(777i)r(m 2 ) 



and the norm factor has $ = ( mi m2 ) Inserting into Eq. (41) we find 

(mi+m2) 

Z 6 (2,s) = - V (-2) mi+ma [mi - m 2 J e Al (m 3 +m 3 ) _ As ( mi+m2) ^ 
' 4 ^— ' mi!m 2 ! mi + m 2 

An important remark is in order. As already discussed in the previous section for the 
case n = 4, one would generally expect that different rapidities cannot have the same 
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values and so one would conclude that m\ and must be of different parities. Indeed 



this is the condition which we used in the Appendix C for the derivation of (73) by 



the direct method. However, as discussed there, this restriction would lead to missing 
states. In fact the states with two strings with ni\ ^ rri2 and with the same parity 
do exist because they are the infinite L limit of deformed strings that at finite L have 
different rapidities. We will discover it more easily using the wedge initial condition 



in Section 5.3 and recover exactly (74), see Eq.(lOl), where the sum is indeed over all 



pairs (mi, 7712). We have not attempted a direct derivation of the contribution of these 
deformed strings. 

It is possible to show that starting from the above expression for Z(2, s), using the 
Airy trick, summing over rrtj and taking the limit t — Y 00 one recovers the second order 
in the expansion of Since it is shown for general n s in the following, we refer to 

section [7] for this calculation. 

This concludes the study of the direct method using fixed size L and the Bethe 
equations. We now turn to the wedge initial condition which, after a few non trivial 
steps, will allow to nicely circumvent the difficulties of the direct method. 



5. The wedge initial condition in half-space 

In this section we report the calculation of the generating function for the wedge initial 
condition. For the symmetric wedge we have only partial results which are given in 



Appendix D[ Conversely, the half-space wedge allows for a full solution that we report 
in the following. 



5.1. Spatial integral of the Bethe eigenf unctions 



Let us consider the (left) half-space wedge model defined by the partition sum (17) with 



w > 0. To compute its moments (39) and their generating function (41) we need the 



spatial integral (38) over the Bethe states. By performing explicit Integration in Eq. 



(|24|) over the negative half-space we find 



= III 



n (1+ 



ic 



n>£>k>l 



Ar, - Ap, 



G 



Ai,..A n 



OO<Xl<X2<..<X„<0 



L jw + i\i + .. + iXj 



where the dependence on w only occurs as an (imaginary) shift of the rapidities. 

It turns out that the sum over permutation can be performed explicitly and can 
be expressed, for any n and any set of rapidities, as a product over pairs of rapidities, 
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a manifestation of the remarkable properties of the BA 

= nl TT i + \ a + \/3 — 2iw 

lla=lV a J 1 < Q</3 < n a p 

In this formula and below we set c = 1. This formula has been also checked with 
mathematica up to large n. 

If we now inject the string solution \j a with a = 1, ..rrij as in Eq. (29), the above 
can be rewritten as a product of intra-string contributions and inter-string pairwise 
contributions in the form 

/W n s 
* M = n!(-2rJ]SSU II V'kM^, (76) 
A — 1 1 ^ A ^ A S », 



i=l l<i<j<n. 



where 



Di , m k = n n '+ a ^+ a ^-. 2 ' w , (77) 



Ko<mj K5<m 



™ _ 1 TT \ + \q + A *,6 ~ 2iw (7R) 

miM (-«)"" n^CV-H^m, V + ^-Mti; ' 1 j 

The intra-string factor . fci can be evaluated for a single string (for convenience we 
choose here the string parameterization Aj )0 = — |(mj + 1 — 2a) + kj, obtained from 
( 29 ) with the replacement a — > rrij + 1 — a) 

s» = ^_fr^- n !{a+b+ 1] - 1 ft 1 fr !{2b) (79) 

-* 11 /(2a) J^ m /(„ + ») 11 /( 2„) 11 /(!+(,) I ' 

1 ^ 1 

with /(x) = —i(m + 1 — x) + 2k — 2iw. Hence we obtain 

gt = ni---H = (zigw = (-ir — , 2=2tt+2ra . (81) 

m,fe r(l-z) T(^ + m) z(2 + l)..(.z + m-l) v 7 

The inter-string factor fc ma fca can be computed similarly 

n» n n /(a + o + i) _jT f(i + b + mi ) 

"mtto.mto- 11 11 /( a + 6) "11 /(1 + 6) 1 ' 

l<a<mi l<b<m 2 J v 7 6=1 ^ v 7 

_ jpg. -igga - 6) + fci + A; 2 - 2zw 

6=1 z 

with /(x) = — i( mi ~^ m2 + 1 — x) + fci + /c2 — 2no. This leads to two equivalent forms 

r (i - z - nn±m) r (i _ z + m±!B2) 
r (i - ^ + 2^2) r (i - z - mf^) 

-ir J, ^-J r !t^J > * = +2*,. (84) 



J-)W _ _V 2 / V 2 / /oo\ 

- Ly mi,fci,m 2 ,fc2 p A , m 1 -m,2 \ p /i mr-mg \ V od / 



r(i- z + ^)r(z+ ., 
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This factor is in general complicated, even for w — 0, although in that case we found 
that its norm simplifies to 

4(fci + k 2 f + (mi + m 2 f 



\D 



io=0 



D 



w=0 



D 



w=0 



mi,/ci,m 2 ,/c 2 I J ^m 1 ,k 1 ,m 2 ,k2 m 1 -k 1 ,m 2 -k 2 A(k\ + k 2 ) 2 + (ffll — m 2 ) 2 

Similarly one also has 



m— 1 



to=0 1 2 r<ii>=0 riii>=0 

— ^m,k ^m,—k 



i qio=i 
I ^m,k 



TT — 

11 zLi-2 _ 



(7=0 



4fc 2 + q 1 



5) 



(86) 



Eq. (|76|) together with Eqs. (|81|) and (|83|) give the complete expression for the 

las 

57) 



space integral of the Bethe states in the half-space model. Note that at large w one has 



1=1 



(-2w) n ' 



A« -» 1, 



id j 



which is the correct normalization of the sharp edge limit since e wy — » 5(y)/w. 



5.2. Generating function 

We have now all ingredients to compute the generating function g(s) 



j _|_^oo Z(n s ,s) 



/7i*=l 



defined in (40) where, from (41) and using (76), we have 



Z(n s , s) = 

oo n s 

e n 

mi,...m„ s =l j=l 



2 m i 



ni 



j Jk 



n 6 



in 

m.i,ki,mj ,kj ' 



where S 1 ™ and Z) w are given in (81) and (83) and we defined 



D w 



(88) 



^9) 



(90) 



where the factor coming from the norm is 

__ 4(fcj - fcj) 2 + (m, - mj) 2 

The main difficulty comes from the double products which preclude independent 
summations over rrii,ki, and can be seen as an interaction between strings. In the 
w — > oo limit (sharp wedge) it simplifies since D w — > 1 and the resulting term can 
be written as a determinant. This allows to decouple the interaction writing the 
determinant as a sum over permutations. For arbitrary w it does not seem that such 
a rewriting exists. In the limit of w — >■ however a determinantal structure (pfaffian) 
also arises, as we will now show. 
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5. 3. Small w limit: study of the poles 

Upon inspection we find that the limit of interest, w = + , is dominated by poles near 
zero momentum in the ki integrations. These are the analog of the states contributing 
in the large but fixed finite size L at w = discussed in Section |4j Those appeared as 
delta functions forcing a set of combinations of sums of kj to vanish for the contributing 
states. Here instead they are poles in the same combinations of kj. The empirical rule 
which appeared there, i.e. that strings appear to pair up in the states with non zero 
contribution, will appear much more clearly here. 

Let us first examine the behavior of the factors S w and D w when w — > + . From 



Eq. (81) there is a divergence at ki = 0, which takes the form 

( -pm,: 

Yim^Aki + 2w) 



qw _ v / i_ „w /qi \ 



where s w is the "regular part" which remains finite for any momentum when w — > 0, 



and is 0(1) when ki ~ w. For the inter-string factor it is clear from Eqs. (84) and (85) 
that when m\ ^ m 2 D^ i k ^ m . k , remains finite for all ki,k 2 as w — > 0. However, for 
mi = m>2 there is a divergence when a pair of momenta vanish (i.e. ki + kj = 0), which 



from Eq. (84) is 



Again we denote d w the "regular part" which remains finite for any momentum when 



w — > 0, and is 0(1) when ki + kj ~ w. By contrast the factor $ in (90) never diverges. 
However it does play an important role: since in each pair of strings a factor $ occurs 
with each D w we should put them together in evaluation of the divergence, i.e. 

We now argue that the only regions in the integral f k k which contribute to 
Z(n s , s) in the limit of small w and of large negative x are obtained by splitting the 
n s = 2N + M strings into iV pairs of strings of (nearly) opposite momenta and M single 
strings of (nearly) zero momentum. More precisely these regions are 

K ~ w, ■■■ h M ~ w , Ki+i + h M +2 ~ w , ■■■ k i M+2 N-i + Kwn ~ w ■ ( 94 ) 

To understand how this works, and the respective role of the limits w — > + and 
x — > — oo let us examine the product 

II s " ' ( 95 ) 

i i<j 

where here we denote for simplicity S { = S^., Dfj = D™. >k . )m . )k . and $^ = <$> m ^ mj)k] . 
We now expand each factor Si = Sf ole + Si and Djj = D p ° le + dij in pole plus regular 
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part according to Eqs. ( 91p2 ). In each of the 2 Ua ^ ns+1 ^ 2 terms of this expansion, each 
string momentum fcj can thus (i) appear in exactly one pole (ii) appear in more than 
one pole (overlapping poles) (iii) appear only in regular parts (unsaturated poles). We 
now argue that: 

- the limit w — > + rules out (ii). 

- the limit x — > — oo rules out (iii). 



Hence the only configurations which contribute in these limits are of the form (94), 
i.e they are such that any given string momentum ki, I = 1, ..n s , appears in exactly one 
pole. The overlapping poles are found to be subdominant for w — > because of phase 
space considerations, and the unsaturated poles decay to zero exponentially as x — > oo. 

Consider first overlapping poles, and focus on the Jacobian in kj integration. 
Let us first check that (94) indeed yields 0(1) momentum integral. One can rescale 
each ki a — > wki a , a = 1,...M and similarly each pair momentum as {ki a _ 1 ,ki a ) — > 
{-K + f K, K + fifeij, a = M + 2, M + 4, ..M + 2JV. This leads to a Jacobian w M+N 
which is compensated by the same factor l/w M+N coming from the poles. 

Suppose now that one of the momenta, say ki, appears both in the pole of a Si 
factor and in one pair factor, say D^. One must thus rescale both ki and kj by w, and 
then one sees from (93 ) that an extra w 2 factor arises due to the $ factor. Even if all the 
3 poles SiSjDij are considered, the jacobian is w 2 while the poles give w~ 3 w 2 = w~ x . 
Hence this contribution is negligible as w — > [§] The same mechanism is at play for 
instance to exclude simultaneous poles in D^DjiDu, since the simultaneous conditions 
ki + kj = 0(w), ki + ki = 0(w), kj + k^ = 0{w) imply ki ~ kj ~ kg = 0(w), hence once 
again the $ (norm) factors in (93) make this contribution subdominant. One can extend 
these arguments to a larger number of overlapping poles with similar conclusions. 

This still does not show that only the configurations (94) remain. To prove this 
we need to show that each string momentum should appear in at least one pole. This 
is the role of the x — > — oo limit. The idea is simple. Let us first integrate over all 
momenta which appear in poles. The remaining integral is 0(w°) hence survives as 
w — > 0, however it now contains momenta which appear in only regular terms. The 



factors exp(— Y^,j m jkjt) can be removed using diffusion kernels as shown in Appendix 
|E|Thus one can first examine the t = integral. Thanks to the factor exp(— i ^ . 
this integral can be performed by contour integration in the sector where all \mkj > 0. 
Examination of poles in this integral shows that it decays exponentially as x — > — oo, 
and this property is not changed by the convolution with diffusion kernels. The general 
argument is then rather intuitive, but giving a general proof is difficult, and for this 
reason we simply examine some examples in |Appendix E| 

Finally, a last property emerges which simplifies even further the calculation. For 



§ this assumes convergence of the rescaled integrals, but examination shows at worse logarithmic 
divergences which do not change the conclusion of being subdominant 
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each pole since we can close the contour in the upper half-plane, we can use the Cauchy 
residue formula which is equivalent to I 



lim — >2n5(k). (96) 

w->o+ ik + w 



This replacement is equivalent, for configurations of the type (94), to perform first the 
w rescaling and then integrate. Replacing each pole by a 5-function, one then obtains 
a x-independent result, which can be shown to equal the limit x — >■ — oo, i.e. the flat 



initial condition for KPZ. This is also discussed in Appendix E 

To summarize, the recipe is as follows: expand the product Yli Si Yli<j D%j in poles 
plus regular part according to 

( — 1 ) mi 

sZm -+ ^y 27 ^) + 4a . (97) 

DZ iM>mj>kj (-l) m< m^ mijm . 2KS(ki + kj) + <C uki)m . )k . , (98) 
with the constraint that each momentum kt appears only in exactly one pole, and 
perform the remaining integrals. The result is the sum of the residues associated to 
configurations where the n s strings split into N pairs of strings of opposite momenta 
with same particle number m and M single strings of zero momentum with all 
distinct number of particles. Of course there are many such partitions, hence many 
contributions, since M and iV can vary from M — 1, ..n s and iV = 0, ..Int[n s /2] with 
the constraint M + 2N = n s . Our claim is that this gives the same result as the finite 
size method of Section [4] as it is directly checked below for all results available from the 
direct method. 

5.4- Examples of configurations 

To illustrate the method, let us consider the simplest class of terms (i) M = n s i.e. n s 
strings with zero momenta (ii) 2N = n s all strings paired. 

5.4-1- All unpaired configurations. The residue of the pole where all ki vanish (pole in 
all Si) reads 

Z(n s ,s)\ N=0 = (99) 

nrrij I 1 \f 

r „(m?-m,')4-Am,s V L ! 



m ! ,0,m,-,0 v 0>m.>0,™i J 



miy£m 2 ^...m„ 3 =l j=l 3 \ Ji l<i<j<n s 

where the summation is over all different m.j. It contains an additional factor | for each 
integration J Note the absence of x dependence. One now uses that for mi 7^ ni2 

D _ ( U min(mum 2 ) ™1 + ™2 $ _ (Wj - m 2 ) 2 

u ' ' 11712 — 171x1 (mi + mi) 2 

|| note that this corrects a misprint in [3S] 
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1 00 1 

*("..«) = as £ (- 2 ) E ' m 'Il^ 

mi,...m n s =l -> 



m^—. — XmjS 

Jo J 



TT f -^\mm(m i ,m i ) \ m i m j\ 
l<j<j'<n s J 

(101) 

where we extended the sum to all rrii, which is possible because of the {mi — rrij) 2 factors 
in the norm. 



5.4-2. All paired string configurations. The other extreme case is for all the strings 
are paired, i.e. 2N = n s , M = 0, assuming n s even. We can label their momenta 
ki, —ki, k 2 , —k 2 , ..k p , —k p , .. with p = 1, ..N and denote their sizes as m 1; m 1; m 2 , m 2 , ... 
We find 

N . „ N 

lJV = l p=l V ""?JJ=" 



oo N N 

mi,...m N =l p=l P Jk Pp=l 



,2(m^— m p ) 8m p fc|A 3 +2Am p s 



JV 
p=l 

X 



n a 



//> /;<> , D° , 



i<p<g<A r 



From Eqs. (86) and (90) we have 



m j 

m,k J m-k^m,k,m-k = Y\_ 4^.2 _|_ ^2 ' 



5° ,52. 



We also see from Eqs. (85) and (90) that 



7)° D c 



mi,fci,m2,fe2 mi,— fei,rri2,— fc2 (f) ' 

^mi.fei.mj ,- fc2 



• (102) 



(103) 



(104) 



and so 



7~)0 r)0 r)0 rjO = <3> <J) flflK 1 ) 

^mp,kp,m q ,k q ^mp,-kp,rn q ,k q ±y m p ,kp,m q ,-k q ly rnp,-kp,rn q ,-k q ^m p ,kp,m q ,k q ^m p ,kp,m q ,-k q j l- LU ' J y' 



which leads to 

1 oo 

ZK = 2iV, s ) = — £ 2 2 ^' 



mi,...mjv=l 



n 



/_iVm p '"P 1 3 

I X / TT 1 c 2(ml-m v )\ — &m p k 2 p\ 3 +2\mpS 

m v J- \ Ak 2 + O 2 

p=l " K P p qr=l P ^ 

n '™ 



A(k p - k q ) 2 + (m p - m q ) 2 A(k p + k q ) 2 + (m p - m q ) 2 



i<p<q<N 



A{k p - k q ) 2 + (m p + m q ) 2 A(k p + k q ) 2 + (m p + m q ) 2 ' 



(106) 
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5. 5. Full generating function 

We now display the complete result for the generating function. We label the 
configurations as follows. There are N strings, p = 1,..N, (nj., ..n p , ..tin) which are 
paired, and M strings which are singles denoted j = 1, ..M with mi, ..rrij, ,.mu- The 
summation over the number of strings can be replaced by independent summation over 
N and M as follows 

N,M>0 y ' 

The partition sum at fixed N, M is now 

oo oo _ M 1 

Z(N,M)= J2 J2 ^ M {-2)^f=^ +2 ^p=^ Yl — - eH-m^^-Xm^ ( 10g ) 

n 1 ,...n N =lm 1 ,...m M =l j=l 

I _ I N r i N np i N 3 

X TT ( l^minfmi.m,) \ m i m 3\ TT / TT TT TT 2(^-n„)^-8n„^A^2An„ S 

J-J- 1 ' m^ + m, J- J- 7 fc n p J- J- J- J- 4fc2 + g 2 11 e 

l<i<j<M 1 3 p=l Jk P Pp=lq=l P y p=l 

t-t 4(fc p - fc 9 ) 2 + (n p - n q f 4(fc p + k q ) 2 + (w p - n g ) 2 

x iJ<:<n ^ - k ^ 2 + k + n «) 2 4 ^ + + k + 

v-iv -r-r 4fc 2 + (n v — ITLj) 2 

x _i Ef=in P (2JV_i)!! rr ^ p 

l<p<JV,l<j<M P v p j; 

where the last line is made of the factors D^^^D^ _ kp>m . >0 = $n p ,fc p , m „o from a 



partial cancellation with the norm factors due to (104). The combinatorial factor 



(2iV — 1)!! = is the number of ways to pair 2N strings. 

This formula allows to recover the moments Z n by expansion in powers of 
exp(— As). We checked it as (i) against the results for the full space obtained in Section 
[4] up to n = 4 (ii) at t = where Z n = 1 is recovered up to n = 12. 

6. Generating function as a Fredholm PfafRan 



We can now take the formula ( 107|108 ) that we just derived as the starting point for 



g(s). We show in this section that it has a determinantal structure, more precisely a 
Pfaffian structure. Hence we first recall the definition of the Pfaffian of a matrix: for 
any antisymmetric matrix A of size 2n x 2n 

n 

aeS 2n ,a(2j-l)<a(2j) 1=1 

with (pf A) 2 = detA. For n = 1 it is pf A = A i2 , for n = 2, pfA = A 12 A U - A 13 A 24: + 
A14A23 and so on, the number of terms of the sum being (2n — 1)!!. 
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6.1. n s even 

Let us first concentrate on the case with an even number of string n s which requires 



in Eq. (108) M to be an even integer and N arbitrary. First we use Schur's Pfaffian 



identity for X u ..X 2n [51] 



Xi + XJ . ...... .V, • ,v, 

J l<i<j<2n 3 



which, for 2n = n s = M + 2N, allows to rewrite the product in Eq. ( 107[ ) as 



n " L i ~ " L j TT 
rn • -I- m • J- J- 



rrii - rrij -p-r 4(k p - k q ) 2 + (n p - n q ) 2 4(k p + k q ) 2 + (n p - n q ) 



{ M nn + mj 4(A; P - k q f + (n p + n q ) 2 A(k p + /c g ) 2 + (n p + n q ) 2 

i<p< 

with 



Ak 2 + (n v + mi) 2 \ 11 2ik v \ X + X, , 

l<p<iV,i<i<M p 1 V p 1 ; \r> I P/ V ' / 7 /; 



X P -i =n p + 2^ p , p = l,...JV, (112) 
Xp =n p -2ik p , p = l,...N, (113) 
Xp+27V = m p , p=l,...M. (114) 

We will need another important Pfaffian identity (valid for M even), checked up 
to large order using mathematica, and which, to our knowledge is new 

J] (-l) min(m! ' mj) sgnK - m 3 ) = pf [(-ln^'^sgnK - m 3 )} MxM , (115) 

l<i<j<M 

valid for any set of integers rrii (we define sgn(0) = 0). 

It is now convenient to treat n p and m p on the same footing. To this aim we define 
a new set of n s integers gj, j = 1, ..n s such that 

SV-i =np, p=l,...N, (116) 
92 P =n p , p=l,...N, (117) 
g p+ 2N = m p , p=l,...M, (118) 

and the n s momentum integration variables Qj, j — 1, ..n s 

Q 2p _! = fcp, p=l,...JV, (119) 
Q 2p = -k p , p = l,...N, (120) 
Q P+ 2N = 0, p=l,...M. (121) 

We want now to sum independently over all the integers gj and integrate over all 
momenta Qj. In order to do so, since g 3 - and Qj are not independent variables, 
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we introduce several delta functions. Putting together the various pieces above and 
inserting in Eq. ( 107[ ) we obtain the intermediate expression 



(2JV-1)!! 



Sl,-fln s =l 
N 



9j 



- 2 ) Ej 93 n / n 2^r^ n ^-^-^-^ 

j=l J Qi q=l 2lC *i + q j=l 



X 



;-l) Ej, = inp pf{TT5jPf|TTTjn 5 



2tt 



j ^77^7^(Q2p + Q-z P -\) \\ 2ir5(Q 2 N+j] 



P =i 



2iQ 2p 



(122) 



3=1 



where pfp^ and pf pTT] > are the Pfaffian in Eqs. (115) and (111) respectively, with 
Xj = g 3 - + iQj. 

Now we can use the Pfaffian identity pf(AtAj-Ay) 2nx2n = (IIi=i A^)pf-4 and rewrite 



2"V 1 TT5|n 27r,5 ^ +J )=pf 



3=1 



4 



_1 WinCffi.Sj) 



sgn(ft - gj) . (123) 



MxM 



Next we can use the permutation symmetry of the integration measure fXj—i Iq 
which allows to make the substitution 



N 



2tt 



(2N - JJ^ 2pi92p i _5(Q 2p + Q 2p _ l} 



-►pf 



27T 

2^ 



5(Q l + Qj] 



(124) 



2Nx2N 



in the integrand since the remainder of the integrand is also a symmetric function in 
the permutations of strings within the 2N paired strings. Note that the matrix inside 
the Pfaffian is indeed antisymmetric thanks to the 5-function (it can also be written 
equivalently as an explicit antisymmetric matrix). 

Now we can further use the full symmetry of the integration measure 
Y[ 2 j=i M ~ n3 Z~2 gj Iq with respect to all n s variables so that again, inside the integrand 
one can replace the product of these two Pfaffian (123) and (124) by a single one 
1 



(2N)\M\ 



Pfp3)Pfp4) 



(125) 



2tt 



, . n m + Qj){ -l)** fttft + ^*(Q0 WiX-l^^Sgate - 9i)) (126) 



using again the symmetry of the remainder of the integrand by permutations of all the 
n s strings. 

Putting all together, we can now write the generating function for n s even in a 
Pfaffian form much more symmetric with respect to all n s strings 



X 

9(s) = 1 + £ —i Z ( n » a ) 



(127) 



n»=l 
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ZK*) = £(-2)^/112^ 

9l >l j=l J Qjq=l ^J^* 



3 

e 



\ — L % / Tig X 71s 



V2i<5j + ft + 2iQj + gj 



:i28i 



This formula has also been checked explicitly to reproduce the same expansion as Eqs. 



(107 108) for n s = 2,4 and up to replica number n = 10. 



6.2. n s odd and general n s 



Now we first need the generalization to M odd of Eqs. (Ill) and (115). In order to 



do this, it is useful to think of the case n s odd as a sum over an addition dummy 
variable g ns +i = tum+i = and over an additional integration over Q ns +i = 0. One 



can then check that Eq. (Ill) remains valid for M odd where the r.h.s. is now a 



(n s + 1) x (n s + 1) Pfafnan with X ns+ i = 0. Equivalently, Eq. (Ill) remains valid for 
M odd if one substitutes there 



,, B V 
Pf(m}^pf _y T 



B^ Wi + gi-nQj-g^ v , = 1 (129) 



13 2iQi + gi + 2iQj + g s ' 

n s +l,n s + l 

in terms of the original two sets of n s variables gj and Qj, i.e. adding a constant unit 
column vector in the Pfafnan. 



A similar reasoning shows that for M odd Eq. (115) must be modified into 



Yl (-l) min(mi ' mj) sgn(m i - m> 



pf 



l<i<j<M 



(-l) min ( m *' m ^sgn(mj -rrij) 1 
-1 T 



M+lxAf+l 

(130) 



where again V = 1 represents the vector with all elements equal to 1. 

The argument of the previous section can then be repeated, introducing delta- 
functions and symmetrizing the integration over the odd n s original variables gj,Qj, 



j = 1, ..n s . By inspection one sees that (|126|) still holds provided one substitutes 

2tt 




U; 



s(Qi 



(131) 



Ai 



2tt 



(n s +l)x{n s +l) 

(2t^ 2 



13 ~ 2iQ^ 1 ^ ™-*> ^ ^ - 9i)- (132) 

We can now transform, for any n s , the product of two above Pfaffians into a single 
Pfafnan using the properties of the Pfaffian of a tensorial product (for n s even) or a 
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slight modification thereof (for n s odd). It reads as follows: consider a 2n s by 2n s 
antisymmetric matrix of the form 

C =(t D T D B \ Aj = UiVj, (133) 

where A and B are two antisymmetric n s by n s matrices and U and V are n s vectors, 
then: 

Pf(<^)2n a ,2n s = pf (-A)n„n,pf ( B )n s ,n s , ™s even, (134) 
pf(<^) 2 n a ,2n s = pf ( _ [/T ^ J pf ( f yT ^ J ,71, odd. (135) 

Applying it to the matrices A, B and vectors U, V defined above we obtain the final 
expression valid both for n s odd and even 

n s r m j _ 9 3 

Z(n s ,s)= J2U II 2/ /, ' "'"''^ ' W 

m,i>lj=l ^ k iq=l 1 3 ^ 1 

( itA ki + fyX- 1 )"* + ^<J(A;i)^(%)(-l) inin(m<,mj) sgii(mi - mj) |(2tt) < 5(A ; 

-1(2^)5(^0 



x pf 



which was previously displayed in Here we restored standard variables Qj — > kj 
and gj — >• to indicate string momenta and number of particle in the string 
respectively. We recall that here and below J k = f^°° y z . This formula has been 
checked against intermediate formula by expanding up to n s = 3,4 and shown to 



reproduce (108) up to n = 10. 



Until now we have not performed the usual shift in the free energy, but we perform 



it from now on. It was discussed in Section (4.2.1) to which we refer for details. It 



amounts to define (62) and focus from now on the generating function associated to 



the variable / = — £ t in Eq. (J3]). 
6.3. Airyzation 

Using the Airy trick (|42|) and shifting y,j — )■ yj + s we get 



Z(n s ,s)= £(-2) E ' m 'II / Ai(y j + S )ll 

m„>l j=\ Jk i^i q=l 



Xmjyj—irrijk^X 3 



2ikj + q 

q= i J 



3 



X pf 



tk 



-5{h + %)(-l) m ^ m! , mj + i?8(ki)8(kj){-\)^™™)s&{mi - mA 7r5(h 

—ir5(kj) 



2n s X 2n s 



3) 2ik i +m i +2ik j +m j / 2 n s x2n s 



where here and below f y = f^dyj. To proceed we will make extensive use of the 
following identity 

XiAijXj X^vA ( A ij U t Vj\ 

Pf 1 -XjUjVi B i3 ) ~ Al " A - pf I -U jVi BlJ J ■ (137) 
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The same identity exists with the role of (A, U) and (B, V) interchanged. It can be 
generalized by formally replacing Aj — > J dviOiivi) or, as below, with summations, e.g. 
Aj — y Ylm a i{ m i)i leading for instance to (for 2n s x 2n s Pfaffians) 



pf 



As 



u iJ ' dvjdjiv^Vjivj) 



-Uj J dVidityAViivi) J dVidVja^Vija^v^B^Vi.Vj) 

Aij 



:i38i 



i=l 



] [ / dviai(vi)]pi 



-UjVi(vi) Bij(vi,Vj) 



i.e. it allows to enter integrations inside the Pfaffian. The same identity again exists 
with the role of (A, U) and (B, V) interchanged. 

There are several ways to proceed from here. In the following we propose one that 
is particularly suited to extract the infinite time limit, but there may be others to arrive 
at different but equivalent forms for the kernels. 

We start by rescaling kj —> kj/X (using Eq. (137)) then shifting jjj — » jjj + 4fc| 
obtaining 



Z(n s ,s)= £(-2)^»*n / M ^ + s + 4 ^ 2 ) n^7A- + q 

f ^5{h + kj){-l) m ^ mim + 7r 2 <5(A; i )<5(fc i )(-l) nlin (' Tl »^)sgn(m i - mj ) 



x pf 



Now we use 
liki + Xnii — 2ikj — Xirij 
2iki + Xirii + 2ikj + Am,- 



-ir5(kj 



7rS(ki) \ 

2iki+\mi—2ikj—Xm,j I 
2ik i +\m i +2ik j +\m j J 2 n s x2n s 



dvidv^Vi - Vj )(d Vj - d v .)e-< 2iki+Xmi) - v ^ +Xm ^ 



3 Jvi>0,Vj>0 



dvidvjS'(vi — Vj)e 



-Vi(2iki+\mi)—Vj (2ikj+\mj) 



(139) 



Vi>0,Vj>0 



since the boundary terms cancel. We can now use Eq. (138) with at(vi) = e M 2 * fc <+ Am <) 
and obtain 

Z(n„ s) = J2 II / M ^ + S + V J + Ak])(-2)^e- 2 ^ +Xm ^ J] — x 



Pf 



gS(h + kj)(-l)^S mi , mj + 7r 2 5(fc i )5(A; J )(-l) min ( m -^)sgn(m i - rnj) nSfcMvj) 

—-Kd(kj)8(vi) 26'(v i -v j ) 



2n s X 2n s 



(140) 



where we have further shifted yj — > yj + Vj. Entering the sums inside the Pfaffian this 
leads to 



Z(n s , s) = Yl Ai (Vj + s + vj + 4/c|)e~ 

1 = 1 -0 
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xpff KniyiihiViW^kjiVj) Kv2{yi,ki,Vi\yj,kj,Vj) \ 

V -K 12 (yj,kj,Vj;yi,ki,Vi) K 2 2(yi,k i ,v i ;y j ,k j ,v j ) J 2n s x2n s ' 

with (we keep the arguments of the kernels implicit) 

Kn = ^{k + k 3 )f ki/x [Ae x ^ + y^} + n^k^k^F^, 2e x ^], (142) 

Ku = 7c5(k t )(e- 2eXy > -l)6{ Vj ), (143) 
K 22 = 25'(v i -v j ), (144) 

where the summations lead to the definition of two functions 

f u ■ yn 1 , z) m_ iF 2 (l;2-2ik,2ik + 2;-z) 

m=l q=l 

F(zi, Zj )= V (-l) min (" i ' m ') S gn(m-m / )^ T ^-(-^r(-^r' (146) 
^— ' ml ml 

m,m'>l 

= sinh(z 2 - Zi) + e~ 22 - e~ 21 

+ / duJ (2y/ z 1 z 2 u)[zi sinh(zi(l — u)) — z 2 sinh(,2 2 (l — u))]. 
Jo 

While the resummation of / is immediate the one for F has been performed from a 



Borel transform as explained in the Appendix F 



The advantage of these manipulations is that the integrations over kj and yj can 
be brought inside the Pfaffian, using again (138), which leads to our main result [^j] 

*<«.,.) =n/ P f( K fr Vi \ ft' V] \) . ("7) 

_,i„,>0 V -Ki 2 {VfVi) K 22 {v i ;V j ) ) 2n s x2n s 

with 

/Jh, -2i{vi-Vj)k 
—dy x dy 2 Ai{y x + s + Vi + 4k 2 )At(y 2 + s + Vj + 4k 2 ) — / fc/A [4e A ^ + ^] 

+ J dy 1 dy2^F[2e Xy \2e Xy2 }At(y 1 + s + Vi )Ai{y 2 + s + Vj ), (148) 

K X2 = K{ Vi )5{v 3 ), K(v t ) = l -j dyAi(y + s + v t )(e~ 2eXy - 1), (149) 
K 22 = 28'(v i -v j ). (150) 

Note that our definition of the Pfaffian of a matrix Kernel K a b(vi,v 2 ) = K(a,vi, b,v 2 ) 
assumes the order 1, v±; 1, v 2 ; ..1, v ns ; 2, v\\ 2, v 2 ; ..2; v ns . Here and below we will use 
indifferently the symbol K(vi,v 2 ) = K n (vi,v 2 ). In evaluations below one should not 
forget that since Ki 2 (vi,v 2 ) is of the form U(vx)V(v 2 ) it cancels in the final result for 
n s even, and it appears linearly for odd n s . In addition the result depends on Kn and 
i^22 only as a function of the product kernel KnK 22 , see Appendix G| 

% correcting a minor misprint in |35j 
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To summarize we have shown that the generating function (|40|) can be written as 



K 12 (Vj;Vi) K 22 (Vi]Vj) /2n s x2n a ' 



n s =l j=l J Vj>0 

Before relating this to a Fredholm Pfaffian let us explicitly evaluate the simplest cases. 



6.4- Evaluation for n s = 1, 2 

Let us first show the one-string contribution, i.e. n s 



Z(l,s) 



K 12 (v,v)=TrK 12 = K(0) 



v>0 



dy 



dv(e~ 2eXv - l)Ai(y + s), (152) 



the evaluation of the Pfaffian for n s = 2 gives simply (see Appendix G ) 
Z(2,s) 



which is identical to the one we obtained in the analysis at fixed n in Section 4.2.1 As 
shown there, this gives the leading asymptotics of g(s) for large s > 0. 

Because K\ 2 is a rank one operator, hence K\ 2 = (TrKi 2 )K 12 and TrK± 2 = K(0) p , 

(153) 



(154) 
sin(2(wj — 

(155) 
(156) 



-TvK ll K 22 = -2 K l0 (v,v) 

Jv>0 

where we have defined the kernel K\o, important for what follows, as 
K w (vi,v 2 ) = d Vl K 11 (v 1 ,v 2 ). 



Notice that since fk is even in k, in ( 148 ) one can replace ^e" 



-2i(vi—Vj)k y 1 



2k 



Vj)k) hence 
Z(2,s) = 



dk 
2^ 

x 



d yi dy 2 K Ai ( yi + s + 4k 2 , y 2 + s + 4k 2 ) 
fk/x[4e x{yi+y2) ] + 2n5(k)^F[2e Xyi ,2e Xy2 



in terms of the Airy Kernel 

K A i(yi, V2) = / Ai{y x + v)Ai(y 2 + v) . 

Jv>0 



(157) 



The product operators read 

(K n K 2 2)(v 1 ,v 2 ) = 2K lQ {v 2: v l ) - 2K n ( Vl ,0)5(v2) , (158) 
(#22#n)K^) = 2K 10 (v u v 2 ) + 25{v 1 )K 11 {^v 2 ) , (159) 

as can be checked by applying to a test function. Hence K 22 K U is equal to 2K W , 
up to a rank one operator which furthermore is nilpotent (i.e. of square zero) since 
#ii (0,0) = 0. 
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6.5. Fredholm Pfaffian 

We now recall the definition of a Fredholm Pfaffian, which we denote as 

I 



<?i(s)=Pf[J + K] 



pf[K(^,U,-)] 2 „ s ,2n s 



Vi>0 



-I 



(160) 



=0 j=l 

where here K is the antisymmetric 2 by 2 matrix kernel of components K a t = K a b(vi,Vj) 
defined in (148 150) and ( 145|146 ) and / is the identity acting in functional space. 
Fredholm Pfaffians are generalizations of Fredholm determinants, for more details see 
e.g. [52llp. 

In ( |l60 ) the symbol pf recalls that the Pfaffian is defined with the order 
1, v\] 2, Vi; 1, v 2 ; 2, v 2 ; -.1, v ns ; 2, v Us . However in our formula (151) we use the other 
ordering 1, v\\ 1, v 2 ; ..1, v ns ; 2, v\\ 2, v 2 ; ..2; v ns , which differs by a permutation a of 2n s 



elements with signature (— 1) 



Pf =(-!)■ 
one has 



"s("s — 1) 



The two definitions are thus related through 



Hence g(s) and gi(s) are closely related, but not identical, i.e. 
g(s) = J2^Z(n s ,s), (161) 

n s s ' 



9i[s 



Now, one useful property of the Fredholm Pfaffian is that it can be written as the square 
root of a Fredholm determinant (FD), namely one has 

I + K t _ K22 
Kn I + K 12 

which is the generalization of the following relation which can be checked for arbitrary 
antisymmetric matrices A, B and arbitrary matrix C: 



gi {sf = Pf [J + K] 2 = Det[I - JK] 



Det 



(163) 



Pf 



.4 



-I-C 



T 



I - 

B 



C 




(164) 



Some basic properties of FD are recalled in Appendix G 



Now we notice that, although not true in general, due to the fact that K\ 2 is a 
rank one operator there also exist a Pfaffian representation for g(s). Going back to the 
way we constructed the Pfaffian form for Z(n s , s) in the previous section we note that 
(i) if n s is even then Z(n s ,s) ~ (KuK 2 2) n ^ 2 while (ii) if n s is odd, due to the fact 



that K\ 2 is a projector, then Z(n s , s) ~ {K\\K 22 
K 22 — > —K 22 provides exactly the factor (— l) ns (™ s 

g(s) = gi{s) K22 ^- K22 . 



|(«.-l)/2tf 12 . 

1 ^ 2 , i.e. from formula (161) 



Hence the formal change 



(165) 



+ equivalently the transformation can be written as ~pi(A) = pi(B T AB) = det(i?)pf (A) where the 
matrix Bij = 5i^ a (j) is a permutation matrix of 2n s elements with dct(£>) = (— 1) s< i ~ . 
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Thus g(s) can be written also as the square root of a Fredholm determinant (FD) 

I + Kf 2 



Det 



(166) 



K 22 

K u I + K 12 

in terms of the Kernels defined in (148 150) and ( 145|146 ). Hence it is also a Fredholm 
Pfaffian, although we will from now on use only the relation to the FD. 

We can perform some further manipulations. Consider the general identity for 
matrices: 



det 



B, 



Sij ~\~ UjV j ■ l - J ; i 

5ij + ViUj 



Aij 



det 



(1 + H(1-A5)-V)) 2 ,(167) 



Bij 
A - A- ■ 

where A and B are arbitrary antisymmetric matrices, and extend it to FD. From 



now on we will be using the standard bra-ket notations. From (166) using that 
-^12(^1,^2) = K(v\)8(v2) is a rank one operator we find 



Det 




(l + {5\(l-K 11 K 22 y'\K)) 



This can be rewritten as 

g(s) = y/Det[I - KuKm](1 + (5\(1 - KuK^K}) 
= v/Deip - K U K 22 }(1 + Tr(l - K lx K 22 y x K X2 ) 
= y/Det[I - K 22 K U }(1 + - KnK n )- x \V). 



In the Appendix H| using (158) we show that 

Det [J - #11^22] = Det[J - K 22 K U ] = Det[l - 2K 10 ], 
and we recall that we denote K w (vi,v 2 ) = d vl K(vi,v 2 ) where K 



K n . 



Appendix H| we also show that 

(#1(1 - KnKu)- 1 ^ = (K\{1 - 2K W )- 1 \5). 
Hence we arrive to the main result for the generating function g(s) 

g(s) = VDet[l-2Ar 10 ](l + (K\(l - 2K 10 )- 1 \5)). 



(168) 

(169) 
(170) 
(171) 

(172) 
In the 

(173) 
(174) 



where we recall the definition of the Kernel Ki in (154) and of the function K in (149). 



7. Large time limit 

We now show that the free energy probability distribution converges to the GOE Tracy 
Widom distribution at large time (i.e long polymer length). In the large time limit (i.e. 
large A 

1 

~2 



(c 2 t/4) 1//3 ) one already sees from (149) that 



lim K[v) 

A— >-oo 



Ai(y + s + v). 



(175) 



y>0 
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Hence from (152) the one-string (n s = 1) partition sum reads 

Ai(2y + s) 



lim Z(l, s) 

A— >oo 



-TVS, 



(176) 



y>0 



where B s (x,y) is the Kernel 

B B (x,y) = e(x)Ai(x + y + s)0(y). (177) 

As shown by Ferrari and Spohn jlS] this is precisely the GOE kernel which enters the 
Fredholm determinant expression for the GOE Tracy Widom distribution 



F^s) = det[J - B s ) = J2 -fin, s) 



n=0 



We will now show that it extends to all n, i.e. that 



lim Z(n s , s) 

A— >+oo 



z(n s ,s) 



-1)™ S / det[S - (a; i ,a; i )]n,xn. 

Jx 1 ,..x„ s 



and equivalently 



lim g(s) = lim Prob(/ > — s) = Fi(s) 

A— >+oo A— >+oo 



(178) 



(179) 



;iso) 



i.e. the free energy probability distribution converges to the GOE Tracy Widom 
distribution. 

Considerable simplifications occur in the expressions fll~48|l49|150p for the FP 
Kernel. First, we can rewrite (175) in the bra-ket notation: 



lim K(v) = — (1\B, 



:i8D 



A— >oo s 2 

involving the GOE Kernel, where (l|(fi) = 1 is a constant vector. Second, we notice 
that the functions fk and F which appear in the finite time Kernel simplify considerably 
in the large A limit 



lim / fc/A (4e Ay ) 



A— >+oo 



lim F(2e x y\2e^)] = 0( yi + y 2 )(%i)0(-2/ 2 ) - B{y 2 )e{-y x )). 

A— >+oo 



;i82) 
;i83) 



The first formula is immediate from ( |145 ) and the second is shown in the Appendix I 



Third, another important property shown in the Appendix is that for A — > +oo: 



2Kiq(vi;v 2 ) 



Ai(y + v\ + s)Ai(y + v 2 + s) - -Ai(s + v{) / Ai(y + s + v 2 ) 

^ Jy>0 



= K Ai (s + vx,s + v 2 ) - Ai(s + Vx)K(v 2 ), 
i.e. it is exactly the Airy Kernel plus a projector and can be written as 

2K-, - d2 



■10 



Bi--\B S 8){IB S 



;i84) 
;i85) 
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We need two identities which are known for matrices but we assume extend readily 
to operators. The so-called matrix determinant lemma 

Det(A+ \u){v\) = (1 + (v|A _1 |u))DeL4, (186) 

for any invertible operator A, and the so-called Sherman-Morrison formula 

(187) 



A- 1 


u)(v 


A- 1 


1 + 


(v 


A' 


-l 


u) 



We use these identities for A = 1 — B 2 S and \u)(v \ = MB S 8)(1B S \. 



Let us now examine the expression (174) for g(s) in the large A limit. From (185) 

Bl 



and using (186) we obtain 
Det(J - 2K W ) = Det 



1 



I — Bg + -\B S 5){1B 



Det 



1 



l-B* 



Now the formula ( 187 ) gives 



{I-2K W )- 1 



B 



2\-l 



1 (1 - B 



2\-l 



B a S)(lB 8 \(l-B_ 



2\-l 



i + l(i| T 



which allows to write the second factor in (174) as 

1 ^ 



1 + (if |(1 - 2K 10 y 1 \5) = 1 - -(l\B s (l - 2K 



;i89) 



1 



■(II 



B x 



1 



(llrffcl^) 



1+B 



2(1 + 



- 1 - 

We can now put all together and write 
= Det(J - B 



gf 



2 = Det[l - 2ATi ](l + (#1(1 - 2K } ^ : 



-10 J 



1+B, 



Det(/-^) T 



1 + I<l|l3^> 

(l + Klli+kl^}) 2 



(190) 
(191) 

(192) 



, + i(l|T J B7l^) + |(l|TTB;l^) : 
where we have used that (1\S) = 1. Now we will use the result of Ferrari and Spohn in 
Ref. 08]: 

1 



Det(J - fi s 
Det(/ + fi s 



1 



-|<5>. 



(193) 



1 + 5, 

Examination of their proof shows that this property holds for any kernel B s (x,y) 
function of x + y only, hence it should also hold for —B s (x + y) which implies that 
one also has 

Det{I + B s ) 1 



Det(J - B s 



(1| 



1 - 5, 



-|<5>. 



(194) 
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Using these two identities we can now simplify ( 192 ) into 

\2 



g(sY = Det(J - B s 



(195) 



which, upon taking the square-root, proves the desired result (180) and the convergence 
to GOE. 



8. Conclusions 



To summarize we have obtained, for arbitrary time, the generating function g(s) for the 
distribution of the free energy of the DP with one free end, equivalently of the height 
of the continuum KPZ interface with a flat initial condition. It takes the form of a 
Fredholm Pfafhan, i.e the square root of a Fredholm determinant. At large time we 
have shown that the distribution of free energy crosses over to the GOE Tracy Widom 
distribution Fi(s). These results should be of importance and testable in experiments. 
It is rewarding that they have been obtained using the replica method, which thus has 
allowed to solve the three main classes (droplet, flat, stationary). 

The method of solution that we found is somewhat indirect, involving a half-space 
wedge model and taking several limits. We have described why this is so, but it would 
be very interesting to find a more direct method to derive the result. It is not clear 
whether this is a purely technical, or more fundamental difficulty On the other hand 
the finite time Kernel that we have found is not very simple. There too, it would be 
interesting to explore whether simpler equivalent expressions could be found. We also 
note that the (half-space and symmetric) wedge model is of interest in itself and deserve 
further studies. Further properties of the solution at finite time, including extracting 
P(f) and numerics are of high interest and left for future studies [4Tj . 
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Appendix A. Statistical tilt symmetry (STS) 

It is easy to show on the replicated Hamiltonian that the moments 

Z n ({xi},t;{yi}, 0) := Z(x 1 t\y 1 ,0)..Z(x n t\y n ,0) , (A.l) 
satisfy the STS symmetry property 

Z n ({ Xi + a + b}, t; { Vi + a}, 0) = e^-* 6 ^-^^}, t; { Vi }, 0) . (A.2) 
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Choosing a = 0, b = x and X{ = STS implies 

Z n ({x}, t; M, 0) = e"^ 2+ ^ yi Z n ({0}, t; {y t }, 0) , (A.3) 
and so 

/ d yi e w ^Z n ({x},t;{ yi },0) = e-^ f e^ w+ ^^Z n ({0} : t; {y t }, 0). (A.4) 
Appendix B. Direct approach n = 4 

Let us give the sum of all contributions, with all integrals displayed explicitly, and then 
discuss each of them separately. We find 

Zi = 8e 5cH - 8e 2c2t + 24c 3 / — , * , c~' kH + cH - 48 / f - - — e' 2 ^^ 

J k (c 2 + 4k 2 ) (c 2 + k 2 ) J k (9c 2 + 4A; 2 ) 

i 48 f c2 0i ~ k lY -2(fef+fe 3 2 )t 

Lm (° 2 + 4k i)( c2 + 4A; i)(c 4 + (ft - k!) 2 + 2c 2 (k 2 + k 2 )) 
where everywhere we denote momenta integrals f k = f f^- Upon explicit momenta 
integration it yields the final result given in the text. We now discuss each term as they 
appear and for each give the spatial integral J ^ and the norm $. 

(i) First term n s — 1 and mi = 4 (one 4-string of zero momentum): 
128 

M/ M = — L , $ = 1. (B.l) 

(ii) Second term, discussed below. 

(iii) Third term n s = 2 and mi = m 2 = m = 2 (2 2-strings, (k + ic/2, k — ic/2, —k + 
ic/2,-k-ic/2)): 

384c k 2 

and we have used m 2 m 2 Y2ki k 2 ~^ ^A m ^ Ik anc ^ a sy mm etry factor of unity. 

(iv) Fourth term n s = 3 (2 1-strings and one 2-string (ki+i^, k\— i|, fc 2 , — 2k± — k 2 ). 
The BE imply k\ = 0. This leads to 

96(9c 2 + 4A; 2 ) 4A; 2 (c 2 + 4fc 2 ) 

fc 2 (c 2 + 4£; 2 ) ' (9c 2 + 4A; 2 ) 2 ' 1 ' 

and there is a factor of 3 from the 3 choices (mi, m 2 , m 3 ) = (2, 1, 1), (1, 2, 1), (1, 1, 2). 

(v) n s = 4, four 1-strings. One checks that BE leave only the possibility k 2 = — ki 
and k 4 = — k 2 which yields 

24c 2 (c 4 + (A; 2 - A; 2 ) 2 + 2c 2 (A; 2 + A; 2 )) r2 .„ . 

= k 2 k 2 (k 2 - k 2 ) 2 L ' (B ' 4) 

while we have 

* = mgkm - klY 

' ~~ (c 2 + 4A; 2 )(c 2 + (h - h) 2 ) 2 (c 2 + 4A|)(c 2 + (h + k 3 ) 2 ) 2 ' { ' ' 



I 



I 
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with a symmetry factor of 3 for three ways to pair the momenta. 

We now come back to the second term. We found this term by considering the 
wedge w > on the infinite space L = +oo in the limit w — > + . It correspond to n s = 2 
and mi = 3, m 2 = 1. As discussed in the text, it is hard to guess from consideration of 
fixed L. In the wedge calculation it corresponds to the contribution 
74, _ f 64 W (8w 2 + 6w + 1) ((fc! - fc 2 ) 2 + 1) 

lmissin9 ~ k M «*i - w + 4) {ki + + (w+ m 



X 



e -4A»(3fc?+fcg-2) (w ^ + ^2 _ 2) + 3^ + 4w 3 _ 5w 2 



(4^ 2 + (2w + l) 2 ) (A: 2 2 + w 2 ) {(ki + k 2 ) 2 + (2w + l) 2 ) ' 
Performing the change of variable ki — > wki which amounts to focus on the pole 
k\ = k-2 = one finds at small w 



Z A \ m is S in g = -32e 8A3 / 1 + O(w) = -8e 8A3 + O(w) . (B.7) 

Appendix C. n s = 2 for direct method 

Consider first 2 m-strings {m\ = m 2 = m, n = 2m) with ki = —k 2 = k. To compute 



J ^ we use formula (49) for the Laplace transform (LT) of the B n (X) factors. After 
inspection we find that we must extract the 1/s 2 pole of the LT, other poles gives 
subdominant contributions after using the BE. In particular all l/s p poles with p > 2, 
i.e. 2 or more zeroes in the set \i p (associated to the set of {A^}), have a zero Ap 
factor. Since all ^ = the expression simplifies into 

n— 1 _. 

* = L £ A P \{— ¥ . (C.l) 

P,Ap^0 i=l ^ 

This allows to compute more easily j *Sf with mathematica up to n = 8 and guess the 
general form given in the text. This form is later established more systematically from 
the wedge initial condition. 

Consider now different strings mi ^ m 2 . We have checked that k\ = k 2 = leading 
to j ~ L 2 is the only possibility, the others give subdominant contributions in L after 
using BE. Here we will consider m x and m 2 to be of opposite parity, otherwise Xj = 
occurs twice which is excluded. The other possibilities are discussed in the text. We 



must thus look for the 1/s 3 in the Laplace transform (49), i.e. two of the must vanish. 

It is then easy to see that only two configurations contribute (i.e. have a non-zero Ap) 
(i) string mi on the left with decreasing imaginary parts and string m 2 on the right 
with decreasing imaginary parts (ii) the reverse. We thus consider 

tc 

Xi = -(mi + 1 - 2i) , l<z<mi, (C.2) 
^m 1+ j = ^-(m 2 + l-2j) , l<j<m 2 , (C.3) 
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which has a Ap = A factor 

A = m 1 \m 2 \ TT TT(1 + 2 .^ 9 . )=m 1 !m 2 ! | mi+m2 | (-ir"(™), 

j-j. J.J. m _ m2 _ i % _)_ 2i \mi — m2 

i=l,?7ii j=l,m,2 

where the and 1712! factors are the intra-string (or single string) contributions using 
the identity ni<i<j<n(H" J=?) = n - Note that this formula makes sense only for strings 
of opposite parity (it diverges otherwise). 

After inverse LT the final result is J ^ ~ L 2 = n\AL 2 /C with 

ni2 m2 mi m\ 

c = ( _ z)mi+m2 - 2 -q ^2 Ami+p yiy, x p = (-o mi+m2 "^(^i)^(^ 2 ), (c.4) 

j=2 p=j i=2 p=i 

m m 

h(m) = JJ ^ |(m + 1 - 2p) = {-l) m i m+ H l - m c- l+m V{m) 2 . (C.5) 

i=2 p=i 

Putting all together it gives the result displayed in the text, using c = — c. 
Appendix D. The wedge initial condition in the full space 



From section 13.11 we want to calculate 

z w ( X ,tr = j2 1 2 - e " / %> (D-i) 



with 



/w poo . n . 

%•■= / (n^j e ^ E ^ iii,iivi/ ^^'-^)- ( d - 2 ) 



J 

We introduce the auxiliary integrals 



£ n , P (Ai,..A n ) = / dxi ^ dXne ix 1 x 1+ ..+i\ n x n+w {x 1 +..+x p )- w {x p+1 +..+x n ) ^ p_gj 

J — OO<Xl<X2<--<Xp<0<X p+ l<..X n < + CO 

and 

B„(Ai, ..A n ) = ^ 5 niP (Ai, ..A n ) , (D.4) 

p=0 

so that 

[d n y* n (y) = n\J2B n (*Pi,-*pJ U (1 - t x )■ (D.5) 

We have 

-Bn,p (Ai, -.A n ) = G p (Ai, ..Ap)G n _ p (— A n , .. — A p+ i), (D.6) 
G p (Ai,..Ap) = / dx!..cix n e (w+lAl)a;i+ - +(w+iAp) ^ 

•/ — oo<xi<X2<..<a;p<0 

1 



(w + i\ l )(2w + iXi + i\ 2 )..(pw + i\i + .. + i\, p ) 
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It should be already evident that the main technical difficulty in this approach is that 
the symmetry of the problem introduces double pieces in the calculations that instead 
will be absent in the half-space model. However, it is a useful model to recover the 
results for the uniform initial condition by simply taking from the start L = +00, and 
considering w — > + at the end. Like this we have obtained the moments Z n , and up 
to n = 4 we have recovered the results of the direct approach. As mentioned there, it 
allowed to identify and obtain the contributions of the deformed strings. We extended 
it also to larger values of n. Since these calculations are not really illuminating, we only 
present here calculations at fixed n s . We present n s = 1 in details and show how the 
various limits are recovered. We will not show details of n s = 2 but we did perform it 



and checked that it reproduces the results displayed in Section 4.2.2 



Appendix D.l. One string partition sum, n s = 1. 

For simplicity we will set c = — c = 1 with no loss of generality (upon rescaling space). 
For a single string the rapidities are A a = — |(m + 1 — 2a) + k and so 

B UtP (w, Ai, ..A„) = G p (w + ik, Ai, ..X p )G n - p (w - ik, -A n , .. - A p+ i) . 

This gives simply 

_ (-2) m Y(l-m-2{w-ik))Y{l-m-2{ik + w)) 

m ' p ~ T(p + l)r(m -p+l)T(l-p-2(w- ik)) T (1 - m + p - 2{ik + w)) ' ^ ' ' 

(-2) m T(l-Aw)T(l-m + 2ik-2w))T(l-m-2ik-2w)) 
m " r(m + l)r (1 - m - Aw) T (1 + 2ik - 2w) Y (1 - 2ik - 2w) ' ^ ' ' 

from which we can write the needed integral 

_ (-2)"T(m + l)r (1 - Aw) T (1 - m + 2ik - 2w) T (1 - m - 2ik - 2w) 
~ r (1 - m - Aw) T (1 + 2ik - 2w) T (1 - 2ik - 2w) • ( • ) 



Inserting the integral in the expressions for Z n we have 
1 1 
n 2 L 



1 1 r w 2 

Z^=^jJ2 ^e~ nkH e^- n) ^ . (D.10) 



k 



Let us first compute for finite w in which case we can replace — > nL A , leading to 

(-2)"r(n)r (1 - 4ii;) (w a_ w)A f kH T{l-2tk-n-2w)T{l + 2ik-n-2w) 
T{l-n-Aw) J k T{l-2ik-2w)T{l + 2ik-2w) 

_ (-2)»r(n)r (1 - Aw) o(n a_ n) « f ^_ nkH ^ 1 



T(l-n-Aw) ° ]3 4fc2 + 0' - 1 + 2w) 2 ' ^ D ' U) 

Using r(— n + x) ~ -, we see that for — > the term j = 1 dominates and acts 
like a, 5(k) function and we get 

2 n 1 2 

z „ = 2 n - l T(n)(n - l)!e (n3 " n) ^ FT ^ = 2"- 1 e (n3 -" ) ^ , (D.12) 

f=2 ~ ^ 



The KPZ equation with flat initial condition and the DP with one free end 45 



recovering Eq. (61) obtained for the uniform initial condition (direct method). 

However this approach allows also to study the wedge problem for arbitrary w and 
not only the flat initial condition. In fact for finite w we have (performing the usual 
shift of the free energy) 

m! 



(—~\) m 

Z(l, s) = J2 ~ j— e" Ams Z^ (D.13) 



m=l 



m=l 



2 m r (1 - Aw) e m s- Xma f _ 4mfc 2 A 3 r (1 - 2ik - m - 2w) T (1 + 2ik - m - 2w) 



mT(l-m- Aw) J k T (1 - 2ik - 2w) V (1 + 2ik - 2w) 



-r (l - Aw) [ Ai(y + Ak 2 + s) 

2 m e mXy p (l - 2i| - m - 2w) V (l + 2l* - 771 - 2w) 



oo 

X 

m=l 



^mT(l-m-Aw) V (l - 2%\ - 2w) Y (l + 2z| - 2w) 

where the rescaling k — >■ k/X has been performed followed by the shift y — >• y + Ak 2 + s. 
Quite surprising the sum can obtained in terms of an hypergeometric function. One 
finds the one-string partition sum 

/dk 

dy—Ai(y + Ak 2 - x)F(X, k, w, y) , (D.14) 

in terms of the function 

F(\,k,w,y) = ~^ /fcX9 1 2 (D.15) 



A (^) 2 + 



w 



> A i 

x e Xy 3 F 3 ({l, 1, 1 + 4^}, {2, 1 - 2i^ + 2w, 1 + 2i\ + 2w}, -2e Xy ). 

A A 

Although the function F seems quite complicated, for fixed A, k and w the large 
y > behavior is very simple 

F(A, fc, y) ~ -y + C + 0(e- iwXy ) + 0( e - eAy ), (D.16) 

as well as F(X, k, w, y) ~ e Xy for large y < 0. 
We can now discuss the various limits: 

(i) very long polymer A — > oo with fixed wedge w. Then we obtain 
Z(l,x)m- dy—yAi{y + Ak 2 + s) = - —w 3/2 Ai{w + s), (D.17) 

since the hypergeometric function behaves as a ln(e Ay ) for large argument and as e Xy 
for small ones. This is exactly the result for the fixed endpoint case (KPZ droplet 
initial condition) which appears as the attractor at large time for fixed w (note that 
the reduced free energy shift in x, \n(2/w)/X is subdominant). 

(ii) flat initial condition w — > at fixed A. Then using that w2+ ^ A p — > Tr8(k/X) 
we find 

F(X, k, w, y) 2n5(k) x ^(e" 26 ^ - 1), (D.18) 
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hence one recovers the result (64) from the direct approach. 

Besides these two limits, the above expression however allows to identify different 
regimes when scaling both w to zero and A to infinity, although this is left for future 
study. 



Appendix E. More on pole structure of half-space model 

Here we illustrate some important properties discussed in the text for the half-space 
model for w — > 0, namely: 

(i) to compute Z(x,t) n at large but finite x,t in the limit w — > one can use the 
substitution — > n5{k) + P^ in the integrals (as described in the text). Neglecting 
terms with finite poles give only exponentially decaying corrections in x. This should 
be useful to describe the universality class Ai 2 — > Ai\, i.e. the flat to wedge class. 

(ii) if furthermore x — > — oo (to describe the flat initial condition class) one can 
use ~~ >" 27r8(k) as discussed in the text. 

We just give some simple illustrative examples with no attempt to construct a 
proof. 



Appendix E.l. No disorder, i.e. n = 1 

It is already instructive to discuss the trivial case without disorder. Equivalently this 
is the case n = 1 in presence of disorder. Consider the half-space model 

Zfot)= f—f dye iky+wy e- k2t - lxk = [ <^^— e -^-^ ( E .l) 
1 ; J 2irJ y<0 y J 2irik + w K J 

< la ze -i f d t_J_ e ^ + a)k^ (K2) 



/Ant J 27r ik + w 
performing the k integration, there is a pole at k — iw. For x + a > the contour can 
be closed in lower half plane and it gives zero. For x + a < the contour can be closed 
in upper half plane and it gives the full residue via the Cauchy formula and so 

Z[x~t)= [ _^ e -^I(i_ sgn ( x + a )) e -(-+ a ) (E.3) 
J V47rt 2 

I x 

~^w^±o+= 2 erfc ^2t I /2' ) ~^ x ^-°° = 1 - 
The two terms correspond to the substitution jA^j = 7r5(k) + Pjfc- The 7r5(k) gives 
the 1/2 while the principal value gives the (odd) sign function. However in the 
limit x — > — oo one sees that for any fixed a the replacement = 2n5(k), which 
corresponds to the Cauchy formula, holds. The conclusion is that for x — > — oo at fixed 
t one can always close the integration contour on the upper half plane, hence use the 
Cauchy formula. In other words the pieces ir8(k) and contribute equally in that 
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limit. Note that to justify this one needs to use convolutions via diffusion kernels. This 
is made more systematic below. 



Appendix E.2. An example with n s = 2 

Here we show on an example that the contributions neglected in the text (case called 



(iii) in the discussion below Eq. (95) indeed decay exponentially. 

Consider n s = 2 and for simplicity the contribution of ni\ = 1 and m 2 = 2 which 

is 

2 2 f 

^(n S ,S)|mi=l,m 2 =2 — 7T / i, l,fc 1 ^2,fc 1 - /J l,fc 1 ,2,fc 2 (I) l,fc 1 ,2,fc 2 e 3 12 3 ) 

Z Jki,k 3 

where we have 

^ = -2^- s '^ - (K4) 

S ^ = 2(ik 2 + wY $2M = ~l + 2ik 2 + 2w' (E ' 5) 
nPO / e n , 3i + 2^! + 2k 2 - tew 

D 1M , 2M - 0, d ltku2M - _ i + 2ki + 2k2 _ 4iw - ( E -6) 

The term S\ ole S 2 ole d\ 2 is the one included in the text. The one left out is 
2 2 



2 



^l,fc 1 S 2,fci«l,fci,2 ) fe*l ) fc 1 ,2,fc2 e 

ki,k 2 

= 2^ e(m 3_ m . ) _|_ Am . s /" 3z + 2fci + 2k 2 - tew . . 

2 " 7^ 2(fei -w)(l + 2^2 + 2^(1 + 2^! + 2i£; 2 + 4w) 1 

4(fei - fe 2 ) 2 + 1 -kft-2k%t-i X (k 1 +2k 2 ) 

4(An - A; 2 ) 2 + 9 

Replacing ^ ^ ^ — >■ |(27r)5(/ci) and then setting to = 0, it simplifies to 

= ^l e K 3 - m j)T2- A "V s / ~0 + 2fc 2 ) -2k\t-ix{2k 2 ) 

"22 " A 2 2{-i + 2k 2 ){-3t + 2ik 2 ) 

= ^i e (m*-m j )±-Xm S 8 f r 1 1 i -2ffi-fa(2fca) /g o\ 

2 2 J k2 [ 2(-i + 2k 2 ) -3z + 2fc 2 J • v • ; 

The integrals at t = are exponentially decaying in x 

1 -2ik 2X = l e -H(i _ sgn(x)), (E.9) 



/ 



2(-z + 2fc 2 



/ * 9 , e" 2 ^ = V 3 N(1- Sg n(x)). (E.10) 
y fc2 -3z + 2A; 2 4 

At finite t they are "broadened" by diffusion but they still decay exponentially: 

j _ e -2klt-ix{2k 2 ) = f da e -£ f J e ~2ik 2 (x+%) 



k2 2{-i + 2k 2 ) J ^/8nt J k2 2(-i + 2fc 2 
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if da _gl r o,_N 1. _,t . / i + A 1. 

= - / ; e 8t e ^ 2 + x ^ = _ ?e x+ 2erfc — - ~ T*e 2 , x ->■ — oo. 

±J a< - 2x V&rt 8 VvVt/ 4 

On the side x — > +oo the decay is even faster, as e~ x ^ 2t \ 
Appendix E.3. A more general argument 

One can make the argument about the limit x — > — oo slightly more general by 
systematic use of diffusion kernel which allows to close the integration contour in the k 
complex plane. We can write 

CO T^s o Til. ■ 

Z(n„ 8 )= Yl ll—e^-^^-^'G^t)^ (E.12) 

mi,...m„ s =lj=l J 

where m = (mi, ..m ns ) and x = (x\,..x ns ) and we have introduced the multi-point 
integrals 

G m (x,f) = n / SZ^e-™^*^ J] (E.13) 

j=l l<i<j<n s 

Writing G m (x, t) as a function of its initial value 

G m (x,t) = T[ / - \ e-^r L ,G m (x + a,t = 0), (E.14) 
J=i J V A7rt / m j 

we have to study simply 

G m (x,0) = n/ S^Ue-*^ II (E.15) 

and in particular its limit when all x^ — > x + aj — > —oo. In that limit this integral can 
be computed by contour integrals closing in Imkj > 0, since all factors S, D are rational 
fractions of the kj and there is an exponential convergence factor e~ lXjTnjkj . We now 
examine the structure of poles and zeros 

i 

S™ poles : kj = iw + -qj, qj = 0, 1, ..m 3 - — 1 , no zeroes, (E.17) 

D™j poles: h + kj = 2iw + iqij, = \\rrii — m,j\, ..X(rrii + m 3 ) — 1. (E.18) 

These are all simple poles and zeroes. They can become multiple however in the 
products, e.g. S^D^Sj has a double pole for = qi + g^. The remaining zeroes 
of Dfj (which do not cancel the poles) correspond to negative = —\\rrii — m,j\ — 
1, .. - \(mi + rrij). 

The poles which have been discussed in the text correspond to qj = and 5^ = 
and their residue lead to the dominant contribution as w — > 0, analyzed in the text. 
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The other poles lead to decaying exponentials as x — > — oo. Note that it is crucial to 
work at fixed t, i.e. the x — > — oo limit should be taken before the large time limit. 



Appendix E.4- Half-space model, n s 



1 



It is instructive to study the one string partition sum n s — 1 of the half-space model at 
fixed w. 



First, using (76) and (81) one has for the single string state with n particles 

An r(l -n-2ik- 2w) _ ^ T(2ik + 2w) 



* = n\ 



(E 19) 

Y{l-2ik-2w) '"" Y{n + 2ik + 2w)' y ' ; 
For large w this gives J w \I/ = (— ) n n! consistent with the result for fixed endpoints. 



Using (39) this leads to 



2 n r(n)e (n3 " n) V 



dfc e _ 4nfc 2 A 3 Y(2ik + 2w) _ inkx 
2ir r(n + 2ik + 2w) 



(E.20) 



with A 3 = t/4 since here we work in units where c = 1. Note that using the diffusion 
propagator of the previous section, one can perform the k integral with poles at 
k = i(w + p/2), p = 0, 1, 2.. and rewrite this as 

da 



Z(x,t) n \ 



n,=l 



^jAixt/n 



-e 4t 



6{-{x + a))e nw{a+x \l + e^ a+x) ) n -\ 



where one recovers (E.3) for n = 1. Now it is clear on this expression that for w — > CP 



and x — > — oo at fixed t one recovers the result from the direct method (61) up to 
exponentially decaying terms. 



Using (E.20) we obtain the partition sum, performing the usual manipulations 



{k — > k/X, Airy trick, shift of y — > y + s + 4/e 2 ), as 

{-l) m 



-\ms 



z(x,ty 



(E.21) 



' ^ dy^Ai(y + 4k 2 + s)^2 



r(2if + 2w) 



XJ "2vr 
Performing the sum we have 

/dk 
dy—Ai(y + 4k 2 + s)F(X,k,w,y) 



\ m r(m + 2%l + 2w) 



e -^j\ (E.22) 



(E.23) 



F(X, k, w, y) = - 2 F 2 ({1, 1}, {2, 1 + 2i- + 2w}, -le^^y 

ik + Xw X 

For large A at fixed w and x one finds F(X, k,w,y) ~ —y0(y) + O(j) hence one 
recovers the fixed endpoint result 

Z(l,s)^- f —yAt(y + 4k 2 -x) = -[ —t 3/2 At(t-x). (E.24) 

Jy>0 27r Jt>0 37r 
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For small w = + and fixed A we can replace l/(ik + Xw) = n5(k) + -4. The delta 
function piece gives 

Z{l,s)\ Spiece = - [ (e- 2eXv - l)Ai(y + s), (E.25) 

4 Jy>0 

independently of x, which is half of the result for the flat space (64). The second piece 
comes from the principal value 

1, s) = -[ dy^ My+ f + s) e^e^ k /\F 2 ({l, 1}, {2, 1 + 2ih, -2<*>e-**'% 
J In ik X 

and can be calculated by antisymetrizing the integral in k. From the above 

considerations it should converge for x — > — oo to the same value ( E.25[ ) and add up to 

the full result (64) as if one had substituted from the start l/{ik + Xw) = 2tx 5{k). This 

statement is equivalent to 

lim - / ^^ e - lxk - tk2 2 F 2 ({l,l},{2,l + 2ik},-2ze- lxk - tk2 ) = ^(e- 2z -l), 
x-^—oo J 2mk Az 

which we have checked via a series expansion. One can see also that at large \x\ the 
piece Zpy(l,s) becomes an antisymmetric function of x. 



jpv 



Appendix F. Summation for the function F 



To compute the function F(zi, z 2 ) defined by the double sum in (146), we note that its 
Borel integral transform satisfies 

Zx Z 2 Z\ — z 2 



dt 1 dt 2 e- tl - t ^F{z 1 t Xl z 2 t 2 ) 



(F.l) 



l 1 + Z\ 1 + Z 2 1 + ZiZ 2 

Denoting s\ = l/z\ and s 2 = l/z 2 , upon the change ti — > U/zi it can be brought in the 
form of a double Laplace transform, hence we have 

1 1 1 s 2 - Si 



F{z 1 ,z 2 )=LT l 



-i 



LT, 



-i 



Sl ^' S2 ^ Z2 s l s 2 l + s 1 l + s 2 l + s l s 2 
si + l)e- 22 + ( Sl - l)s 



LT 



-i 



s\ + l)e 



'i 



8 X (Sl - 1) 



-e Zl ~ Z2 - e 



si 

-*i + e -*z + J (2yfz^z~ 2 ) + / dtJ (2 v ^ 2 l)(e z ^ t - e 



t—Z\ 1 



(F.2) 

(F.3) 
(F.4) 



we have used the (formal) result that the LT of e TZ f Q z dte ±l g{t) is g(s)/(s ± 1). This 
function is antisymmetric in z%, z 2 and to make it more explicit we make some change 



of variable which brings it in the form (146) given in the text. One can then check 



explicitly that it admits the double Taylor expansion given in the text. Note that 
although it looks like divergent at large z±,z 2 there are massive cancellations and the 



function F is in fact bounded and well behaved. This is further discussed in Appendix 

m 
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Appendix G. Fredholm determinants and Pfaffian 

Appendix G.l. Series expansion of a Fredholm determinant. 

It is useful to recall the series expansion of a Fredholm determinant of Kernel M(x, y) 
where / is the identity: 



Det(J + M) — e Tlhl{I+M) = 1 + Tr M + -((Tr M) 2 - Tr M 2 ) 
+ ^ ( (Tr M) 3 - 3 Tr M Tr M 2 + 2 Tr M 3 ) 



(G.l) 



+ — ((Tr My - 6(Tr My Tr M + 3(Tr My + 8 Tr M Tr M A — 6 Tr M ) + . . . . 

which can also be used unambiguously for a matrix operator M a b(x, y) = M(a, x; b, y), 
a = 1,2, b = 1,2, e.g. it gives 



Det(7 + M) = Tr M n + Tr M 22 



(G.2) 



+ -((Tr M u + Tr M 22 ) 2 - (Tr M 2 X + Tr M 2 2 + 2 Tr M 12 M 21 )) + 0(M 3 ). 

Appendix G.2. Generating function g(s) in terms of pfaffians: series expansion. 
Consider now the generating function studied in the text 

oo ^ n s „ 

9( S ) = Yl — Z K) I [ / Pf[K(Ui,Uj)]2n s ,2n, , 

with the definition of the pfaffian of a matrix Kernel K ab (v 1 ,V2) = K(a,Vi,b,v 2 ) with 
the order l,vi, l,v 2 ; ..l,v ns ;2,v 1 ;2,v 2 ] -2;v na , namely 

Kwiv^vi, 
Kuiv^vx) 

/ K 11 (v 1 ,v 2 ) Ki 2 (vi,vi) K 12 (v u v 2 ) \ 

K n (v 2 ,vi) K 12 (v 2 ,vi) K 12 (v 2 ,v 2 ) 

-K 12 (v!,vi) -K^v^v-l) ^22(^1,^2) 

y -K 12 (vi,v 2 ) -K 12 (v 2 ,v 2 ) K 22 (v 2 ,vi} J 

1 + TrK 12 - \{{TrK 12 f - TrK 2 l2 + TrK u K 22 ) + 0(K 3 ) 



Pf[J + K] := 1 + 



Jv!>o y 



+ 



2! 



V!>0,v 2 >0 



(G.3) 



+ 0(K 3 ) 



and, more systematically 
Z(l)=TrK 12 , 

Z{2) = TvK 2 l2 - (TrK 12 ) 2 - Tr K n K 22 , 

Z(3) = -(Tr ^ 12 ) 3 + 3 Tr K 12 Tr fT 2 2 - 2 Tr K\ 2 - 3 Tr X 12 Tr K n K 22 
+ 6Tr K U K 22 K 12 , 



(G.4) 
(G.5) 
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Z{A) = Tr[K 12 } 4 - 6 Tr[^ 12 ] 2 Ti[K 2 2 ] + 3 Ti[K 2 2 ] 2 + 8 Tr[K 12 ] Tr[K* 2 ] - 6 Trfi^] 
+ QTt[K 12 ] 2 Tt[K 11 K 22 } - 6Tr[K 2 2 ] Tr[K n K 22 \ + 3 Ti[K n K 22 } 2 - 6 Ti[(K n K 22 ) 2 } 
- 24Tr[K 12 ] Tt[K u K 22 K 12 \ + 2ATi[K u K 22 K 2 2 ] + 12Tr[K u KT 2 K 22 K 12 \, 

which simplifies drastically when K\ 2 is of rank one, as it is the case in this paper, 
i.e. K\ 2 = Tr{K\ 2 )K\ 2 or Ki 2 (vi,v 2 ) = U{yi)V{y 2 ). Using also the antisymmetry of 
K u , K 12 we have 

Z(l) = TrK 12 , (G.6) 
Z(2) = -TrK n K 22 , 

Z(3) = -3 Tr K 12 Tr ^11^22 + 6 Tr K U K 22 K 12 , 
Z(A) = 3Tt[K u K 22 } 2 - 6Tr[(K n ir 22 ) 2 ], 

which was used in the text. 



Appendix G.2.1. Series expansion of a Fredholm determinant, gi(s). Now consider the 
matrix operator 



M = — JK = " 7 ^" = V 2 r , (G.7) 





and let us consider the square root of the following FD considered in the text, which 
is a Fredholm Pfaffian. We will write its series expansion to second order. Using Eq. 
(G.l ) one finds 

Pf [J + K] = v/Detp - JK] = y/Dct(I + M) (G.8) 

= e 5 Trln ( 1+M ) = 1 + -TrM+ - ((TrM) 2 — 2 Tr M 2 ) + 0(M 3 ) 
2 8 

= 1 + Tr K 12 + ^((TtK 12 ) 2 - TtK 2 2 + Tr K n K l2 ) + 0(K 3 ). (G.9) 



Now we see that this result differs by a sign from the above expression (G.4) for the 
term n s = 2. In general we expect that it will differ by (— i)M n «- 1 )/ 2 . This is because 
the connection from FD to FP relies on the order 1, v%\ 2, V\\ 1, v 2 ; 2, v 2 ; ..1, v n& ; 2, v ns , 
as we now explain. Note that in our case, since (Tri^ 12 ) 2 — Ti K 2 2 = we can simply 
change K 22 — > —K 22 to switch from one formula to the other one, as explained in the 
text. 



Appendix G.3. More pfaffians and Grassmann path integrals 

It is useful to use Grassmann integrals to manipulate pfaffian and determinant, and 
Grassmann path integrals for Fredholm determinants and Fredholm Pfaffian. Consider 
Grassman integral with rjf = 0, J dr] = and J drjrj = 1. We will use convention 
J drji..drj n rji..rj n = 1, the other convention J drji..drj n rj n ..rji = 1 differs by (— l)"-^ 1 )/ 2 
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since r} n ..r]i = (— l) n< ^ n ~ l ^ 2 rji..rj n (the same identities below apply usually with e x — > 
e~ x ). Note that the linear change of variables gives a jacobian = det -1 . 

One has for an antisymmetric matrix A of size n = 2p even, and an arbitrary 
matrix B one has 

j f[dr)ie* mA w = p f[ y4 ] ) J dfj l dr) 1 df) 2 dr) 2 ..dfj n dr) n e fiiBi M = detB. (G.IO) 



This allows to show easily that if A is antisymmetric of even size n = 2p, writing 

Vi = + «A) an d fji = -^(oci- ifa) one has 

det[A] = J df] X dri l dfi 2 dri2.df] n dri n e TllAl3 ' n3 (G.ll) 
da x d(i x da 2 d^ 2 ..da n d^ a ^ a ^ A ^ = pf[A] 2 , 



since the cross term i(ai(3j — ^a^A^ = (ai(3j + aj(3i)Aij = 0. From J dfjidrji = 1 one 
gets that dfjidrji = dai ^ 1 . The sign of the permutation ai(3i..a n /3 n — > aia 2 ..a n /3i(3 2 ..(3 n 
is (— l) n (™ -1 )/ 2 . For n = 2p this is (— l) p = i n . Note that the determinant is indeed zero 
if A has odd dimension. 



Now let us consider the formula (164) given in the text for block matrices made of 
arbitrary antisymmetric matrices A, B and an arbitrary matrix C. It corresponds to 
the other definition of the pfaffian, noted pf. To clarify let us consider the grassman 
integral representation 

det ( I+ f T ~* c j = J \\drji{i)drji{i)dfj 2 {i)drj 2 {i)e s (G.12) 
+ X^^'i^W^O') + C(j,i)fji(i)rji(j) - B(i,j)fji(i)rj 2 (j) + A(i,j)f} 2 (i)r]i(j). 



'j 

It turns out that the transformation 



ViU) l 




< i i 


mU) J 




y-i 1 


ViU) ' 




1 1 i 


m(J) J 







«i(j) -ipi{j) 



(G.14) 
(G.15) 



«i0') + 

«a0') 

decouples the a (z) from the /5 (z). One can check that it leads to the square of the 
pfaffian pf . This is because there are two possible Grassmann measure 

JJ^ (z) = d7 ?1 (l)d^(l)%(2)tf7^(2)...d7 ? i(n a )d^(n s ) (c), (G.16) 
= ^i(l)^i(2)....d»7i(n»)d%(l)d%(2)...d%(n,) (d), (G.17) 
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the first one allows for the above change of variable and corresponds to pf and the 
second to pf . 

Finally a Fredholm determinant can be seen as a one dimensional path integral 
f JJ^(x)rfr/(x)e^^ {x,)(<5{x - y)+ ^ {:c ' 2/)),?(y) = Det[I + K\. (G.18) 

It can be expanded around the free theory, which allows an easy derivation of its 
expansion in terms of partial determinants 



i r 

Y[df](x)dr](x)e^ dx ^ x)ri{x) J2^( I dxdyf](x)K(x,y))r](y)) n 



n=0 



(G.19) 



If 

= / dx 1 ..dx n dy 1 ..dy n K(x 1 ,y 1 )..K(x n ,y n )(f](x 1 )T](y 1 )..f](x n )r](y n )) 

n=0 ' J 
°° 1 f 

= 5^ — / dx 1 ..dx n ^2(-l) a K(x 1 ,x ail) )..K(x n ,y a(n) ) (G.20) 

n=0 ' ^ a 
°° 1 f 

= dx 1 ..dx n det[K(x i ,x j )] nxn . (G.21) 

We have used fj Xl r) a ( Xl )...fj Xp r) a ( Xp ) = (-l) tr fj Xl r] Xl ...fj Xp T] Xp which, up to multiplication by 
Yix(l + VxVx) yields the desired result. 

Similarly one can define from the Grasmann path integral the FD of a matrix 
operator and Fredholm Pfaffians. We will not pursue further but indicate again that 
there are the two possible choices for the Grassmann measure J7 a dr) xa (x) which singles 
out the correct definition of the FP. 

Appendix H. Manipulations on g(s): some details 



Here we show the formulas (172) and (173) in the main text. Let us use the matrix 



determinant Lemma (186) applied to operators, with A = 1 — 2K W , \u) = 2\S) and 

(H.l) 



v) = \K(0, .)). To prove (172) using (158) we need to show 



(v\A- l \u) = 2(K(0, .)|(1 - 2K W )- 1 \5) = 0. 
Since (K(0, .)\S) = K(0, 0) = 0, it is thus sufficient to show that for all n > 1 

Q n := j K(0 :Vl )(K 1Q ) n (vi,0) = 0, 

which we now prove. 
For n = 1 we have 

.1 

t>i>0 z 



(H.2) 



Qi 



(H.3) 
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since K(0, 0) = and K(vi,v 2 ) vanishes when any of the argument goes to +oo. Here 
we denote K^iv^v^) = d l Vi dl 2 K{vi,v 2 ). For n = 2 we have 

J Vl,V2>0 

= - [ K(0,v 1 )K^ 1 \v 1 ,v 2 )K(v 2 ,0) 

J vi,V2>0 

= I K^ 1 \0,v 1 )K^ 1 \v 1 ,v 2 )K(v 2 ,0) = -Q 2 , (H.4) 

Jvi,v 2 >0 

where in the second line we have integrated by parts first over v 2 then over v\. The 
integration on either v 2 or V\ has no boundary term since K(v, 0) = —K(0, v) vanishes 
both at v — and v — oo. In the last line we have used antisymmetry of K and that 
KW(a,b) = -K^(b,a). Hence Q 2 = 0. 

The same manipulations can be pursued for n > 3, although now one has to keep 
track of the boundary terms. One writes 

Qn = f iT(0,t;0^ (1 ' H^1^2)^ (1 ' H^^3)..^ (1 ' H^l,^)^ (1 ' 0) K,0), 
J vi,..v„>0 

and we can now integrate by part all variables starting from v n 

Jvi,..V n >0 

- [K(0, vjKlWfa 1*)]%^ if K^\v u v l+1 )K(v n , 0) 

1=2 

n-1 p-2 



+E(-i)^(o^i)n^ (i ' o) ^,% + i)^ (i ' o) (v^^ (o,i) K 1 %i)] 




p=2 j=l 
n-1 

x J] K^ 0,1 \vi, vi + i)K(v n , 0) 
i=p+i 

n-2 

+ (-l)»Jf(0 )Wl )n ^ (1 ' 0) fo, [^ (1 ' 0) K-i, 0)]^= °° } , (H.5) 

J'=l 

where we have collected all boundary terms. Now using the antisymmetry of K, the 
property K^ 0,1 \a,b) = — K^'°\b, a) and relabeling the integration variables, the above 
sum can be rewritten as (each term in the order they appear) 

n-1 

Qn = -Qn + + ^Q p _iQ n _ p + 0, (H.6) 
p=2 

where the second and last line above vanish due to K(0, 0) = and the third can be 
expressed from Qk with k < n. Hence we have proved by recurrence that Q n = for 
all n > 1. 



The KPZ equation with flat initial condition and the DP with one free end 



56 



To prove (173) we now use the Sherman-Morrison formula (|187|) on the form (|158|) 
(1 - K 22 K lx y x = (1 - 2KX0)- 1 - 2 



[1-2K W )- 1 \5)(K(0,.)\(1-2K 



10 



1 + 2(K(0, .)|(1-2A 10 )- 1 |(5) 



(H.7) 



multiplying left by (K\ and right by \5) we immediately obtain (173) since we just 



proved that the second term vanishes from (H.l). 



Appendix I. Simplification of F[2e Xyi , 2e Xy2 ] in the large A limit 

Here we analyze the large A behaviour of 

h x (y u y 2 ) = F[2e Xyi ,2e Xy2 ] = sinh(2e Ay2 - 2e Xyi ) + e" 2 ^ 2 - e" 2 ^ 1 (LI) 

+ I duJ (Ae^ x{yi+y2) ^u) [2e Xyi sinh(2(l - u)e Xyi ) - 2e Xy2 sinh(2(l - u)e Xy2 )}. (1.2) 
Jo 

We first manipulate this formula to make more apparent the fact that this function is 
bounded. For that one can use the identity 



> i 
duJ (Ae^ +y2 ^)e Xy2 e 2 ^ eXv2 = t^ e -^\ 



to rewrite 



h x ( yi ,y 2 ) = e-^ V2 -e-^ yi 



duJ (4e^ +y2 ^)[e Xy2 e 2 ^ eXy2 - e ^ e ^-)^] 

+ [\uJ (Ae^ +y2 ^)[e Xy2 e- 2{l - u)e " y2 _ e^e" 2 ^^ 1 )]. 
Jo 



Let us write it as 



hxiyuto) = e~ 2eXv2 - e" 2eAH1 + I 2 - h + J 2 - J u 



+oo 



)e Xy i 



J; 



duJ (4e^ x{yi+y2) ^)e Xy 'e 2{1 - u 



duJ (Ae^ +y2 ^)e X ^e- 2 ^ eXVj 



For j = 1, 2. Note that |Jo(<£)| < 1 for z > hence we have the exact bounds 



(I.3) 
(L4) 

(1.5) 

(1.6) 
(1.7) 

(1.8) 
(1.9) 



Hence the function h\ is bounded. 

Let us now examine the limit A — > oo. With no restriction we can set yi < y 2 . 
There are four distinct regions: 

(i) yi < and y 2 < then h 00 (y 1 ,y 2 ) = 0. 

(ii) yi < and y 2 > and yi + y 2 < . 
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For yi + y 2 < the argument of the Bessel function in J± j2 can be set to zero. The 
integral can then be done exactly, which leads to 

J 2 = ±(l- e -**»)«i, J 1 = l(l- e -^)«0 (1.10) 
Next we can rewrite 

r+oo i 

I 2 = dvJ (4:ei x{yi+y2 Wl + ve- x yi)e- 2v « -, (1.11) 
Jo 2 

r+oo 

h= dwJ (4e-^ yi+y2) VTT^).e Xyi e~ 2weXvi « (1.12) 
Jo 

J 2 being clearly is dominated by the argument of J (z) near z = and converges to 
1/2. Evaluation of 7i with mathematica indicates that it vanishes. The final result in 
that region is then 

h 00 (y 1 ,y 2 ) = 0-1 + ^- + ^- = 0. (1.13) 

(iii) y 1 < and y 2 > and y 1 + y 2 > 

In that region one finds J\ ~ 0, ii ~ and one can write 

J 2 = f duJ (Ae^ +y ^Vl^u)e Xy2 e- 2ueXv2 (1.14) 
Jo 

= / rfMJ (4el Afel+ ^ ) V / l - Me- A ^) e -2« ~ J (4e^ 1+ ^) _ ^ 
Jo 2 

and 

/+oo 
duJ Q {Ae^ + ^^)e x ^e 2 ^ eV2 

r+oo 

= / ^Jo(4e5 A ^ 1 +^)^r+^)e A? ' 2 e- 2ueA! ' 2 
Jo 

r+oo i 

= / duJ {Ae-2 Xiyi+y2 Wl + ue-^)e- 2u « - J (4e^ 1+ ^) « 0. (1.15) 

The final result in that region is then 

M!/i,2/2) = 0- 1 + + + + = -1. (1.16) 

(iv) yi > and y 2 > there we can safely find that all integrals are zero hence 

Mj/i»sfc) = °- 

Using the antisymmetry of h\, the final conclusion is that 

hooiyiM) = %i + y2)(0( yi )9(-y 2 ) - 9(y 2 )9(- yi )), (1.17) 
as indicated in the text. 
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Appendix J. Useful Airy identities 

A useful identity involving Airy functions is [55] 

Ai(u)Ai(v) = ^J *^(2 2/3 (£ 2 + U -±^))e^ 

which allows to prove the "duality" identity 
dQ 



■■»)€ 



2tt 



Ai(a + 4Q 2 )Ai(b + 4Q 2 )e 2icQ 



d% A .,a + b + c 



2tt 



Ai( 



,a + b — c 



+ 4?)e 



2\ e 2i{a-b)(, 
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(J.l) 



(J.2) 



Appendix K. Simplifications of the Kernel K\q in the large A limit 



Consider the definition (148) of the Kernel Ku in the large A limit (i.e large time) and 
insert the limit forms ( 182|l83 ). The result can be written as the sum Ku = K u + K n 
with 



Ai(yt + s + Vi + AQ 2 )Ai{y 2 + s + Vj + 4Q 2 



2/1,2/2, 



-2i(vi-Vj 



X 



2iQ 



4 



poo pO 

\ dyi dy 2 Ai(y 1 + s + Vi)Ai(y 2 + s + v j )-(v i *+v j ) 
Jo J- yi 



(K.1) 



Rewriting K n and using the duality identity of Appendix J we obtain 

__ 1 
~ 2 



dc Ai(yi + s + Vi + AQ 2 )Ai(y 2 + s + Vj + 4Q 2 ) 

yi,V2,Q J-(vi-Vj) 
x e- 2icQ 9( yi + y 2 ) 



1 f Vi ~ Vj , [ d£ 
dc 



(Vi-Vj) 



2n 



2/1,2/2 



Ai( yi+y2+v 2 i+v i +c +s +±e) 



x M{ yi+V2+V l + V 3 C + s + ^2 )e 2 l{ v^ ]+yi -y 2 )i e{yi + ^ 



I fVi-Vj 



dc / Ai(y + 



Vi + Vi + c 



+ s)Ai(y + 



Vi + Va — c 



+ s) 



-(vi-vj) Jy>0 

the integration over y\ — y 2 produces a factor tt6(£). Computing the derivative, using 
integration by parts one finds: 

2 Jy>0 



d Vi K u ( Vl - Vj ] 



Ai(y + Vi + s)Ai(y + Vj + s) 



(K.2) 
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1 /" ■ ~~ ' ' 

8 



dcAi{- — j- + s)Ai(- — j- + s) 



-(Vi-Vj) 

and note that (d Vi — d Vj )Ku(vi',Vj) = KAi(vi,Vj). The derivative of the second part of 
K\\ reads 

d Vi Kl 1 (v i ,v j ) = ~ I [Ai(y + s + Vi)Ai(-y + s + Vj) - Ai(-y + s + Vi)Ai(y + s + Vj)} 

4 7j/>0 



— ^-Ai(s + Vi) j Ai(y + s + vj 

4 Jy>0 



(K.3) 



Adding (K.2) and (K.3) one obtains Ki . One notes that then the second line in (K.2) 



exactly cancels the first line in (K.3). To see it we note that the second line in (K.2) 
can be rewritten as 



while the first line in (K.3) can be rewritten as 
1 



4 



daAi(a + s)Ai(vi + Vj — a + s)[9(vi < a < Vj) — 6(vj < a < Vi)], 



(K.4) 



(K.5) 



and the second piece in (K.2) exactly cancels the first line in (K.3). This proves the 



formula (184) and (185) for K w in the large A limit. 



Appendix L. Explicit check for n s = 2,3,4 of convergence to the GOE 



Here we check the property (|179|) explicitly up to n s = 4. The general expansion of a 

(L.l) 



Fredholm determinant is given in (G.l ) with M = —B s which yields 

-TrB s , 



z{l,s) 

z(2,s) = (TrB s ) 2 -TrB 2 s , 

z(3, s) = -(TiB s ) 3 + 3TtB s TtB 2 s - 2Tr5 s 3 , 

z(4, s) = (Trfi s ) 4 - 6(Tr J B,) 2 Tr J B s 2 + 3(Tr5 2 ) 2 + 8T rJ B s Trfi s 3 - 6TrS s 4 . 



On the other hand, calculation of the Pfaffian is given in (G.6) and reads 

Z(l,s) = TrK 12 
Z(2 7 s) = -TrK n K 22 , 

Z{3, s) = -3TtK 12 TtK u K 22 + 6TyK u K 22 K 12 = 3Z(1, s)Z{2, s) + 6TrK n K 22J fC 12 , 
Z(4, s) = 3(TiK 11 K 22 ) 2 - 6Ti(K u K 22 ) 2 = 3Z{2, s) 2 - 6Tr(A- n ^ 22 ) 2 . 

For n s = 2 we have 

Z(2,s) = -2 j Kio(v,v)= j (-Ai(y + v + s) 2 + \ai(s +v)Ai(y + v + s)). 

Jv>0 Jv>0,y>0 ^ 



(L.2) 



(L.3) 
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On the other hand 

Ai(x x + x 2 + s) 2 , (L.4) 



TtB s = / Ai(2x + s), TrB 2 = I 

Jx>0 Jxi 



>0,x 2 >0 



while the first term in (L.3) is clearly —TrB 2 , checking that the second is (TyB s ) 2 
requires some combinatorics. Note that it does not require any property of the Airy 
function, and, as for all values of n s below, replacing Ai(z) by any function would 
give the same result. Hence to check it the simplest method is to write formally 



Ai( Zi ) = J t f(ti)t 



or equivalently to replace each Ai(zi) — > e tiZ \ integrate and 



symmetrize over the % Being true for any set of U it is true for any function. We find 
(Tr5 s ) 2 ->■ e - si - tl+t2) /(4tit 2 ) while the second term in {lJ| yields e ~ s( * 1+ * 2) /(2t 2 (t 1 + t 2 )) 
which under symmetrization is identical. Hence we find Z(2, s) = z(2, s). 



To check n s = 3, from (L.2) we simply need to show that 
6TrK n K 22 K n = 2(TiB s ) 3 - 2TrB 3 s , 



with 
TrB 3 



Ai(xi + x 2 + s)Ai(x 2 + x 3 + s)Ai(x 3 + x\ + s). 



(L.5) 
(L.6) 



xi>0,a;2>0,a;3>0 



We have, from (158) 
6TtK u K 22 K 12 = 6 

= -3 



(K u K 22 ){v 2 ,v 1 )K{v 1 )S{v 2 ) = 12 / K w { Vl ^)k{ Vl ) 



i)i>0,f 2 >0 



i-i >0 



Ai( Vl + s + v 1 )[Ai{y 2 + + s)Ai(y 2 + s) - -Ai(s + v^Aifa + s)}. 

Ui>0,yi>0,s/ 2 >0 

Let us again use the exponential trick replacing each Ai(zi) — > e~ Zlt \ after integration 
we find 

3(t 2 - t 3 ) 



12 / K(v 1 )K 10 (v u 0) = sym tlMi 



vi>0 



t3 2t 1 t 3 (ti + t 2 )(t 2 + t3) 



1 



= 2(TtB s ) 3 -2TtB 3 s 



(L.7) 
(L.8) 



4tlt 2 t 3 (tl+t 2 )(t 2 + t3)(tl+t3) 

upon symmetrization (we are not showing the trivial e _s ^ 1+ * 2+ * 3 ^ factor). Let us note 
that 

(TrB s y = 



2Pt 1 ..U 



Tr B p = . 

S (h + t 2 )(t 2 + t 3 )..(*n-l + t n )(tn + h) 



To check n s = 4, we see from (L.2) that we only need to show that 
- 6 TriKnK^) 2 = -2(Tr B S ) A + 8 Tr B s Tr B 3 - 6 Tr B A S . 



(L.9) 



(L.10) 



We have, from (158|) 
Tr(K n K 22 ) 2 = 



Kio(v 1 ,v 2 )Kio(v 2 ,v 1 ) 



tll>0,«2>0 



K(v,0)K 10 (v,0). 



v>0 
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Note that the second term can be integrated by part and vanishes. 
We also have 

2K 10 (v 1 ,v 2 ) = h ~ h e-* 1 " 1 "* 2 " 2 , (L.ll) 
m 21 2t 2 (t 1 +t 2 ) ' V 1 

hence after multiplying with the same with exchange of v\,v 2 and variables £3^4 and 
integration over Vi,v 2 we find 

ATrK 2 =qvtti {h ~ t 2 ){h - U) 



Multiplying by —6 this should be compared with the r.h.s of ( |L.10 ) 



1 1 

+ 8sym- 



l^hhhu {h + 1 2 ) (t 2 + 1 3 ) {t 3 + u) {u + h) 2t x (t 2 + 1 3 ) (t 3 + U) {u + 1 2 ) ' 

and we find again that it is equal, i.e. Z(A, s) = z(4, s). 
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